arXiv:1506.06902v2 [math-ph] 3 Apr 2017 


A BISPECTRAL g-HYPERGEOMETRIC BASIS FOR 
A CLASS OF QUANTUM INTEGRABLE MODELS 


PASCAL BASEILHAC AND XAVIER MARTIN 


Abstract. For the class of quantum integrable models generated from the q—Onsager algebra, a basis of bispectral 
multivariable q—orthogonal polynomials is exhibited. In a first part, it is shown that the multivariable Askey-Wilson 
polynomials with N variables and Af + 3 parameters introduced by Gasper and Rahman |GR041| generate a family 
of infinite dimensional modules for the q—Onsager algebra, whose fundamental generators are realized in terms of 
the multivariable q—difference and difference operators proposed by Iliev |I108| . Raising and lowering operators 
extending those of Sahi ISaOTI are also constructed. In a second part, finite dimensional modules are constructed 
and studied for a certain class of parameters and if the N variables belong to a discrete support. In this case, 
the bispectral property finds a natural interpretation within the framework of tridiagonal pairs. In a third part, 
eigenfunctions of the q—Dolan-Grady hierarchy are considered in the polynomial basis. In particular, invariant 
subspaces are identified for certain conditions generalizing Nepomechie’s relations. In a fourth part, the analysis is 
extended to the special case q = 1. This framework provides a q—hypergeometric formulation of quantum integrable 
models such as the open XXZ spin chain with generic integrable boundary conditions (q yf 1). 
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1. Introduction 

In the context of mathematical physics, studying the relation between g—hypergeometric special functions in 
one or several variables, related orthogonal polynomials, and the representation theory of quantum algebras is 
an important field of research. From different point of views detailed in the text, this subject is challenging and 
deserves to be further explored for several reasons. 

Starting from the basic example of the quantum integrable harmonic oscillator whose eigenfunctions are the Her- 
mite polynomials, generalizations involving a coupled system of multiple harmonic oscillators have been extensively 
studied in the literature (see e.g. |SJII) 1. In particular, orthogonal polynomials of the Askey scheme |AW85) are 
eigenfunctions of Hamiltonians characterized by a certain tridiagonal interaction structure describing ‘triangular’ 
two-dimensional spin lattices associated with the XY spin chain with nearest neighbor non-homogeneous couplings 
and non-zero magnetic field (see [GMVZI^ and references therein). According to the interaction considered, the 
parameters entering in the definition of the g—hypergeometric eigenfunctions are restricted. For the simplest ex¬ 
amples of models, corresponding eigenfunctions are expressed in terms of one-variable polynomials such as the 
g— Racah, g— Hahn, g— Meixner or g— Krawtchouk ones. Importantly, thanks to the explicit expressions of the 
Hamiltonian’s eigenfunctions in terms of g— hypergeometric functions and using orthogonality properties, closed 
form expressions for correlation functions have been derived in various cases (see e.g. [S.1 1 1] 1. In the last twenty 
years, such models have found applications in different areas, ranging from low dimensional quantum physics and 
statistical mechanics to game theory, also impacting the analysis of strongly correlated systems or Markov chains. 
Indeed, they provide examples of transmission channels for short distance quantum communication in which the 
so-called ‘perfect state transfer’ occurs [CDELO.'i] . In the context of Markov chains, orthogonal polynomials also 
appear in the spectral representation of the transition probability matrix [AGRll] . Applications of Askey-Wilson 
polynomials in the context of solvable quantum mechanics have also been considered. See for instance [OS 14) and 
references therein. 

Multivariable extensions of the basic hypergeometric orthogonal polynomials have also arisen in the literature 
as eigenfunctions of mutually commuting differential or difference operators. In this vein, Macdonald-Koornwinder 
polynomials with N variables and parameters g, t |M881 IM03I IKo921 ISIIl ISa97] are one possible generalization of 
the Askey-Wilson polynomials. They provide a basis of eigenfunctions for relativisitic integrable generalizations of 
the Calogero-Sutherland-Moser system associated with different root systems |R95llD9^ID95llFSII] . Also, certain 
degenerations of the Macdonald-Koornwinder polynomials known as g— Whittaker functions (t = 0), respectively 
Hall-Littlewood polynomials (g = 0), appear as eigenfunctions of the Hamiltonian of the open Toda chain with 
different types of boundary conditions, respectively of the repulsive delta Bose gas |DE1411 IDEI421 [006] . The limit 
A^ —>■ oo of Macdonald polynomials has also been considered recently |NS14] , and for g — >■ 1 Macdonald polynomials 
reduce to Jack polynomials which describe eigenfunctions of the quantum Benjamin-Ono model [NSI3j . In the 
context of conformal field theory and integrable perturbations, explicit formulae for singular vectors of the Virasoro 
[MY95j and g— Virasoro algebra jSSAFR.O^ are given by Jack and Macdonald polynomials, respectively. Let us 
also mention that Macdonald-Koornwinder polynomials provide a basis for a class of solutions of the boundary 
g— Knizhnik-Zamolodchikov equations [SVI3) . 

From the point of view of representation theory and harmonic analysis on quantum algebras, g— hypergeometric 
orthogonal polynomials can be understood as natural extensions of the hypergeometric functions associated with 
Lie groups. For SUq{2), it is known that matrix elements of irreducible representations can be expressed in terms of 
little g— Jacobi polynomials. Askey-Wilson polynomials also found an interpretation as ‘zonal spherical functions’ 
on SUq{2) (for a review, see e.g. |Koe96) l and as overlap coefficients between eigenbasis that diagonalize certain 
elements of the so-called Zhedanov’s algebra [Z91] (the ‘Cartesian’ g— deformation of SI 2 , see also |Fai90[[HPK96j L 
also known as the Askey-Wilson algebra and often denoted AW{i) in the literature^. Macdonald-Koornwinder 
multivariable polynomials of various types also appear as zonal spherical functions in the context of the reflection 
equation algebra |L021 INoDS97j . In the representation theory of affine Hecke algebras, it is also known that each 
family of orthogonal polynomials can be associated with an irreducible affine root system [C051 IC95j . 


^Note that the Askey-Wilson algebra first appeared in the context of Leonard pairs under the name ’Leonard algebra’. See ITer87l 
Definition 3.1]. 
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Besides the Macdonald-Koornwinder generalization of the Askey-Wilson polynomials, a different type of mul¬ 
tivariable generalization of the Askey-Wilson polynomials has been proposecH by Gasper and Rahman [GR041| . 
These polynomials can be seen as g—analogs of the multivariable polynomials previously introduced by Tratnik 
[T89| (see also [GerI07) l. In the present paper, Gasper-Rahman and Tratnik’s multivariable polynomials and their 
connection with the g—Onsager algebra [Ter032l |Bas04] and the theory of quantum integrable systems will be at 
the center of our attention. 

Let us now describe the content of this paper. First, besides the two-variable cas^l, to our knowledge the relation 
between Gasper-Rahman multivariable generalization of Askey-Wilson polynomials |GR041) and the representation 
theory of quantum algebras has not been considered up to now. However, recall that Askey-Wilson polynomials 
with four parameters and g appear in the representation theory of the g— Onsager algebra [Ter032| . generating a 
simple example of infinite dimensional irreducible module. Also, observe that by restricting the four parameters and 
the variable to a discrete support, the Askey-Wilson polynomials generate a finite dimensional irreducible module. 
In this case, the well-known bispectral property of the Askey-Wilson polynomials finds a natural interpretation 
within the theory of Leonard pairs |Ter87[ ITer03I] . In the first part of this paper (Sections 2,3), we will fill a 
gap in the literature: we first show that Gasper-Rahman multivariable polynomials generate a family of infinite 
dimensional module of the g— Onsager algebra. We also construct raising and lowering operators which generalize 
those for the Askey-Wilson polynomials, proposed by [San7) . Secondly, restricting the parameters and choosing a 
discrete support for the N variables, the bispectral property of Gasper-Rahman multivariable polynomials will be 
interpreted within the theory of tridiagonal pairs |Ter032) . Parallel to the case g ^ 1, a similar analysis is handled 
in Section 5 for the special case g = 1. These are the first main results of this paper. 

Now, in the context of quantum integrable systems establishing a bridge between the representation the¬ 
ory of the g— Onsager algebra and Gasper-Rahman polynomials opens new perspectives. Indeed, a systematic 
g— hypergeometric formulation of the eigenfunctions of all mutually commuting quantities of a given integrable 
model in terms of known orthogonal polynomials of one or several variables is highly desirable. In addition to the 
examples mentionned at the beginning of the introduction in which polynomials of the Askey scheme, Macdonald- 
Koornwinder generalizations or limiting cases appear, there are two important families of quantum integrable 
systems for which a g— hypergeometric formulation remains essentially, up to now, unrevealed: 

(i) The first class follows from the Faddeev-Reshetikhin-Takhtadjan and Sklyanin’s (FRTS) framework, where 
the Hamiltonians of quantum integrable models are derived from a generating function of mutually commuting 
quantities called the transfer matrix. The transfer matrix is built from solutions of the Yang-Baxter and reflection 
equations which constitute the basic and main ingredients in the construction (see [Sk88] for details). In this frame¬ 
work, solutions of the Yang-Baxter and reflection equations can be interpreted as intertwiners of finite dimensional 
irreducible representations of quantum algebras. Heisenberg spin chains with periodic, twisted or open integrable 
boundary conditions fall among the well-known examples built in the FRTS framework, associated with Uq{sl 2 ) or 
higher rank generalizations. 

(ii) The second class of quantum integrable systems is generated from the Onsager algebra |044l IDav90| or its 
g— deformation [Bas04) . For instance, the Hamiltonians of the Ising [044] . superintegrable chiral Potts [GeR,85) . 
XY models [Ar90] and generalizations [AS9I) are in one-to-one correspondence with the simplest element of an 
Abelian subalgebra of the Onsager algebra (the so-called Dolan-Grady hierarchy) acting on a finite dimensional 
irreducible module. Similarly, the Hamiltonian of the open XXZ chain with generic boundary condition^ can be 
expressed as a linear combination of mutually commuting quantities which generate the so-called g— Dolan-Grady 
hierarchy [Bas04| , acting on a finite or infinite dimensional irreducible module (in the thermodynamic or continuum 
limit). 

For a large class of quantum integrable systems which have been studied in the framework (i) and/or (ii), it is in 
general not clear how to express the eigenfunctions of the transfer matrix in terms of g— hypergeometric functions 
and related polynomials. However, in recent years it appeared that quantum integrable models with boundaries 
associated with a Uq{sl 2 ) i?—matrix, and AT—matrices solutions of the reflection equation, can be alternatively 
generated from a g— deformed analog of the Onsager algebra [Bas04] or, more generally, one of its higher rank 

^Note that the two-variable case was constructed in IKDJ96I . 

^For q ^ 1, see IKDJ96I . For q = 1, two-variable Racah polynomials of Tratnik IT89I arise in the representation theory of an 
extension of the Racah algebra QR{3), see Em. 

"^For triangular, diagonal or special Uq{sl 2 ) invariant boundary conditions, quotients of the q—Onsager algebra have to be considered. 
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generalization [BB09) . In these cases, the models admit a presentation either within the FRTS framework (i) 
or within the g—Onsager framework (ii). For the simplest Uq{sl 2 ) case, the correspondence between the two 
frameworks has been studied in details |BK05) . 

In Section 4, we will consider a subclass of quantum integrable systems which possess two alternative presen¬ 
tations (i) and (ii) for generic values of q. The Hamiltonian H - or more generally the corresponding transfer 
matrix - can be written as a linear combination of mutually commuting quantities {l 2 fe-i-i|^ = 0,1,2,...} that form 
a g—deformed analog of the Dolan-Grady hierarchy. Namely, 

N 

(1.1) H ='^h^khk+i + K with [\ 2 k+i,hi+i] = 0, k,lGN, 

k^O 

where fc = 0,1, 2,... and Hq are model-dependent^ scalar functions. Let wq, wi, g+, g_ be scalars. The elements 
bfe-i-i are given bj0: 

(1-2) l2fc-i-i = woW_fe -I- wiWfc+i + g+Gfe+i -I- g-Gk+i, 

where W_fc, W^+i, G^+i, G^+i generate the current algebra of the g— Onsager algebra |BSh09] . For all known 
examples of the literature considered up to now, W_fc,Wfc+i (resp. Gfc_|_i,Gfc+i) are polynomials of total degree 
2A: -I-1 (resp. 2k + 2) in the fundamental generators Wq, Wi of the g—Onsager algebra with defining relations shown 
later in Their explicit expressions can be computed recursively using the algorithm given in [BBlOj . For 

instance, the elements Gi,W_i,G 2 read: 

(1.3) Gi=[Wi,Wo]^, 

W_i = - ((g2 -H g-2)WoWiWo - WoWi - WiWg) + Wi , 

^2 = ^ - (g3 + g)Wi2w2 + (g-^ - g3)(WoWi2Wo + WiW^Wi) 

p{q^ + q-^) V 

+ 2g-i)WoWiWoWi + (g^ + g^ + 2g)WiWoWiWo -f p(g - q-^){^l + Wi^)) . 

The expressions for the elements Gi,W 2 ,G 2 are obtained from Gi,W_i,G 2 by exchanging Wq O Wi. Note that 
the vector space on which Wo,Wi act is specified according to the favourite model considered. Quotients of the 
g—Onsager algebra may be considered as well (see e.g. |BB12) L In this ‘g—Onsager framework’, the explicit 
connection between Gasper-Rahman multivariable polynomials and the representation theory of the g—Onsager 
algebra mentionned above gives a straightforward access to a g—hypergeometric formulation of the eigenfunctions 
for this class of models, which is the second main result of this paper. 

The paper is organized as follows. Below, we introduce notations that will be used thoughout this paper. In 
Section 2, the orthogonal system of multivariable polynomials of Gasper and Rahman [GR,041| is recalled. Then, we 
introduce g—difference and difference operators that essentially follow from Iliev’s work |I108j . and are diagonalized 
by Gasper-Rahman polynomials. Studying the structure of the operators, new infinite dimensional modules of the 
g—Onsager algebra are exhibited in terms of these polynomials. In particular, explicit realizations of the standard 
generators of the g—Onsager algebra in terms of multivariable g—difference and difference operators are given. 
Generalizing the construction of Sahi [Sa07] to N^-variables, explicit realizations of raising and lowering operators 
are also given. In Section 3, restricting the {N -|- 3)—parameter space and the variables to a discrete support, finite 
dimensional modules of the g—Onsager algebra are constructed and described. In this case, the corresponding 
difference-difference bispectral problem is naturally interpreted within the theory of tridiagonal pairs. Based on 
these results, the spectral problem for the g—Dolan-Grady hierarchy of mutually commuting quantities (13 is 
reconsidered in Section 4, where the eigenfunctions are shown to admit an expansion in terms of Gasper-Rahman 
polynomials. In particular, subspaces that are invariant under the action of the Abelian subalgebra generated by 
the elements US) are identified, extending the analysis of [BKQ7j . These subspaces are characterized by certain 

^For the Ising or superintegrable Potts model (q = 1), = 5^0? h-Q = 0 with Kronecker symbol 5}zi- For the open XXZ chain, the 

explicit expressions for and are given in |BK07| . 

^The parameters cjq, cji, are model-dependent. For the Ising and superintegrable chiral Potts model {q = 1), cjq = 1, tJi is 

the free magnetic field coupled with the system. For the open XXZ chain, are associated with the boundary parameters 

IBK07I . 
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relations that already appeared in the literature on boundary integrable models, named ‘Nepomechie’s relations’. 
In Section 5, the construction of a hypergeometric basis for the special case g = 1 is handled: infinite and finite 
dimensional modules of the Onsager algebra and a subclass of tridiagonal algebras, respectively, are generated 
by Tratnik’s multivariable polynomials [T89) of Krawtchouk and Racah type, respectively. In this special case 
g = 1, standard generators are expressed in terms of Iliev-Geronimo’s difference operators proposed in |GerI07) . 
Perspectives for the analysis of quantum integrable models and potential generalizations of the present framework 
are briefly described in the last Section. 


1.1. Notations. In this paper, we denote q the deformation parameter, assumed not to be a root of unity. We 
will use the standard g-shifted factorials [GR 042) : 

n—1 oo k 

(a; g)„ = (1 - ag'^), (a; g)oo = (1 - ag^) and (ai, 02, ..., Ofc; g)„ = ]^(aj; g)„. 

fe=o fe=o j=l 


Define Vz = <C[zf ^as the ring of Laurent polynomials in the variables zi, 22 ,..., zat with complex 
coefficients. For each of the variables Zj, j = 1, 2,..., iV, introduce an involution Ij of Vz such that Ij{zj) = zj^ and 
Ij{zk) = Zk for k = 1,2 ,..., N and k ^ j. Then, we will denote Vx = C[xi,X 2 , ■■■, xat ] the subring of Vz consisting 
of polynomials in the variables Xi = {zi + z ~^)/2 with complex coefficients. 

Let {ei, 62 ,..., Cat} be the canonical basis for . We also defincQ the g-shift, forward and backward difference 
operators in the j-ih coordinate acting on functions /(z) = /(zi, Z 2 ,..., Zn), respectively, such that: 

EzJ[z) = f{zi,Z2,...,q^Zj,...,ZN), 

Azjiz) = (Ez^-l)f{z), 

Vzjiz) = il-E-^)fiz). 

Similarly, will denote the forward shift in the variable Uk- For every function /(n) = f(jii,n 2 , ■ ■ ■ ,'UAr) we 
have EriJin) = /(n + 6 /,). 

The standard multi-index notation will be also used. For instance, if u = (ui, V 2 , ■ ■ ■, vn) £ then 


Z = Z-^'-Zo 




e:=eze:i 


■E, 


zq 


= (zig"Gz 2 g"G...,zwg'^"). 


Let v = (ui, U 2 ,..., uat) £ {0, ±1}'^\{0}'^. Following [1108] . we define = max{i'j,0) and = —inin{i'j,0). We 
denote I'^ the composition of the involutions Ij corresponding to the non-zero coordinates of v~. Also, we denote 
W\ = Vl + 1^2 + + UaT. 


2. The g—O nsager algebra and g —difference operators 

Tridiagonal algebras have been introduced and studied in |Ter93[ IITTer99l ITer032) , where they first appeared in 
the context of P— and Q—polynomial association schemes. A tridiagonal algebra is an associative algebra with unit 
which consists of two generators A and A* called the standard generators (see Section 3, Definition 13.21) . In general, 
the defining relations depend on five scalars p, p*, 7 , 7 * and /3 [Ter032) . In this Section, we will focus on the so-called 
reduced parameter sequence 7 = 0, 7 * = 0, /3 = g^ -I- g“^ and p = p* which exhibits all the interesting properties 
and can be extended to more general parameter sequences when g 7 ^ 1. In the literature, the corresponding algebra 
is called the g—Onsager algebra, in view of its close relationship to the Onsager algebra [044] and the Dolan-Grady 
relations |DoG82] . 

Definition 2.1 f [Ter032) IBas04[ lBB09] b The q—Onsager algebra Oq{sl 2 ) is the associative algebra with unit and 
standard generators Wq, Wi subject to the following relation^ 

(2.1) [Wo, [Wo, [Wo, Wi],],-i] = p[Wo, Wi] , [Wi, [Wi, [Wi, Wo],],-i] = p[W,, Wo] . 

^Compared to the notation Eq,Zj in [1108I , one has Ezj = E^2 ^. ■ 

®The g—commutator [J^, y]^ = qXY — q~^YX, where q is the deformation parameter, is introduced. 
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Remark 2.1. For p = 0 the relations i2.1\} reduce to the q—Serre relations of Uq{sl 2 )- For q = 1, p = IQ they 
coincide with the Dolan-Grady relations |Do(;82| . 

In [Ter032] . an infinite dimensional irreducible module of the g—Onsager algebra based on the Askey-Wilson 
polynomials was constructed. On the vector space of all polynomials in the single variable x, it was shown that 
the elements Wq, Wi act, respectively, as a multiplicative factor in the variable x and as the Askey-Wilson second 
order g—difference operator. Alternatively, in the basis of Askey-Wilson polynomials, the elements Wq, Wi act as a 
tridiagonal matrix and diagonal matrix, respectively. Thus, the well-known g—difference equation and three-term 
recurrence relations of the Askey-Wilson polynomials can be interpreted as a bispectral g—difference-difference 
problem for the elements Wq, Wi. Note that similar objects appear in previous works [GLZ92[IGH962[lNoS00| and 
that the connection with the g—Onsager algebra is emphasized in [Ter032j . 

In this Section, we construct a family of infinite dimensional modules of the g—Onsager algebra on the vector 
space Vx based on the Gasper-Rahman polynomials in the multivariables Xi,X 2 , ■■■,xn with coefficients in C, thus 
generalizing the example of |Ter032) . As we will show, to the action of Wo,Wi, one can associate a bispectral 
system of coupled recurrence relations and multivariable g—difference equations that is solved by the polynomials. 
In particular, the case A^ = 2 is described in detail. Starting from the bispectral system of equations, we then 
extend the analysis of Sahi [SaO?] that holds for = 1 to generic values oi N: N distinct pairs of raising/lowering 
operators for the Gasper-Rahman polynomials are obtained, expressed in terms of the standard generators of the 
g—Onsager algebra. 


2.1. Bispectral multivariable Askey-Wilson g— orthogonal polynomials. First, let us recall the definition 
of the multivariable polynomials introduced by Gasper and Rahman in [GR041) using the notations of [1108] . 
As multivariable generalizations of the well-known Askey-Wilson polynomials, these polynomials are obtained by 
applying a Gram-Schmidt process. Let g, a, b, c, d denote nonzero scalars in C. To avoid degenerate situations, we 
assume g is not a root of unity, and that none of ab, ac, ad^ be, bd, cd, abed is an integral power of g. For n = 0,1, 2,... 
the Askey-Wilson polynomials are defined by [AW85] : 


looi f r {ab,ac,ad;q)n 

(2.2) pn{x; a, b, c, d) = --- 403 


b abedq" 

ab, ac, ad 




1. In 

with X = -[z —). 

2 z 


By definition, they satisfy a second-order g—difference equation and a three-term recurrence relation of the form: 

(2.3) p.{z,z~'^))pn{x) = (g“”-I-a6cdg”“^)pn(a:), 

(2.4) bnPn+l{x) anPn{x) CnPn-lix) = XPn{x), 

where piz), p{z, z~^) are rational functions of z, z~^ and an, bn, Cn are rational functions of g" and of the parameters 
a, b, c, d. Their explicit expressions can be found in [AW85] . 

If a,b,c,d are real, or occur in complex conjugate pairs if complex, and are such that max(|a|, |&|, |c|, |(i|) < 1, 
the polynomials satisfy the orthogonality relatioiu: 


[ pn{x;a,b, c,d)pm{x-,a,b,c, dx = 5n 

271 Vl — x^ 


{abedq'^ g)„(a6cdg^"; g) 


' {q'^^^,abq^, acq^, adq^, beq^, bdq^, cdq'^] q)c 


where 


w(z; a, b, c, d) 


_^;g)oo_ 

{az, az~^, bz, bz~^, cz, cz~^,dz, dz~^; g)oo 


In the limit g —t 1 and various limiting processes on the parameters a, b, c, d, one can obtain from the Askey- 
Wilson polynomials all the classical polynomials of continuous and discrete arguments: Wilson, Racah, Hahn, 
Jacobi, Krawtehouk, etc |KS96) . 

In |T89| . a multivariable generalization of the Racah-Wilson polynomials was constructed. From that point of 
view, Gasper-Rahman polynomials |GR041) are g—analogs of the multivariable polynomials introduced by Tratnik. 
They constitute a multivariable generalization of the Askey-Wilson polynomials (12.2L defined as follows: 


®If a > 1, b, c, d are real or one is real and the other two are complex conjugates, max(|b|, |c|, |d|) < 1 and the pairwise products of 
ab, ac, ad, be, bd have modulus less than one, then the polynomials satisfy a different orthogonality relation IKS96I eq. (3.1.3)]. 
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Definition 2.2 (See |I1081 IGR.041) 1. Define Pn{x] a, b,c,d) = Pn{x;a,b,c,d)\q^q 2 from / fg.gl) . Introduce the N -{-i 
nonzero parameters ao,ai, ...,aN +2 and denote Xj = ^{zj + z~^). Identify Zq = ao and z^+i = aN+ 2 - The 
N—variable generalization of the Askey-Wilson polynomials is defined by: 


(2.5) 


N 


gW({n},{x},{a}) = Y[p^.{xj;ajq 

i=i 




^3 „297,_i _ 

, - Zj+i, 


<f3+l-l ^ 
^7 + 1 '' 


where ‘Tlj = ni + n 2 + ... + Uj and ^Tto = 0. 


An analog of the orthogonality property of the Askey-Wilson polynomials can be exhibited as follows. Introduce 
the measure on [— 1 , 1 ]-^: 

( 2 . 6 , 


^ dr 

jj- axj 


and the inner product: 
(2.7) 


if, 9 )= [ f{x)g[x)dp[x). 

Jf-l.ll'V 


*7 


Following [I1081 Theorem 2.1], by induction on N it is straightforward to show the following theorem: 
Theorem 2.1 (See |I108) . Theorem 2.1). Assume: 

(2.8) 0 < loN+il < lawl < • • • < |q;i| < ?7im(l, laop), 


|q:n+ 


j— < |q:a+ 2 | < 


\aN\ 


OlN\ loAT+ll 

The multivariable polynomials ({n}, {a;}, {a}) are orthogonal with respect to the inner product Namely: 

(2.9) (Q(^^({n}, {x}, {a}), g(^)({m}, {a;}, {a})) = 

with 

r-r I 9 N I 9 joo 


= n 

k 


1 


Vl (q2(nfc + l)^ 


For further convenience, normalized multivariable polynomials can be introduced following [1108) . Assume the 
complex variables {z} = {zi, Z 2 , ■■■, zn) G (C*)^, the complex variable {n\ = (ni, n 2 ,..., u-at) G (C mod 
and the parameters {a} = (oo, oi,..., aAr+ 2 ) G (C*)'^’*'^. 

Definition 2.3 (See |I108j l. The normalized N—variable Gasper-Rahman mutlivariable polynomials is defined by: 


( 2 . 10 ) QW(M,{x},{a}) = ,- 


{aN+iaN+2, “"+1?'^+= ; g 2 )^^ 


gW(M,{x},{a}). 


Compared with (12.5F the normalization factor is chosen as follows. Consider the following change of variables 
and parameters: 

aoaAf+iaAf+29 r o Ar , 1 ~ 

ao i-t ao = ao, a^ H> a^ =- tor j = 1 , 2 ,..., A-|- 1 , aAr +2 a 7 v +2 =-, 

o:N-i-2-j Oioq 


Zi I—>• Z-i = 


_ aN-i-2-j ^2<nN+i- 


aoq 


^2n, ^ ^2n, ^ aN+l-jZN+l-j j = l, 2 ,...,A. 

aNAr2-jZN-^2-j 


^^Then, determines uniquely nj. 
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The map f : {n,z,a) i->- {h,z,a) defines an involntion on (C mod x (C*)^ x (C*)^+^ (see [11081 

Lemma 5.1]). Applying Sear’s formula (cf. |GR042j . page 49, eq. (2.10.4)), according to the factorized struc¬ 
ture of the Gasper-Rahman polynomials (12.51) in terms of Askey-Wilson polynomials it is possible to show: 

Theorem 2.2 (See [1108] . Theorem 5.3). The normalized Gasper-Rahman multivariable polynomials are invariant 
under the action of the involution f; 

QW(M,{4,{a}) = QW({n},{i},{a}). 

By analogy with the Askey-Wilson polynomials that solve the bispectral problem (|2.3I) - (I2.4I1 . the normalized 
Gasper-Rahman multivariable polynomials (12.101) also solve a bispectral problem associated with g—difference and 
difference operators. This is the subject of the next subsection. 


2.2. Iliev’s g— difference and difference operators. Following [1108] . we now introduce a commutative algebra 
Az of g—difference operators and a commutative algebra An of difference operators. Let 

P, = C(, Z 2 ,..., \ (BJ',..., ] 

denote the associative algebra of g—difference operators with rational functions of zi, Z 2 , zn as coefficients. The 
commutative subalgebra Az of Vz is generated by N algebraically independent g—difference operators = 

1, 2,..., iV} that are defined as follows: 


Definition 2.4 (See |I108j '). Let v = (i^i, P 2 , l'a) G {0, ±1}'^\{0}'^. Let ..., he the nonzero compo¬ 
nents of V with 1 < < 12 < • • • < is < A" and s > 1. For v € {0,1}'^\{0}^, denote 


^ /I N nfc=2(l 




)( 1 - 


q OLi,Zi,Zi 




Qin 




^ 0!-N+lO^N+2Zi^ aN+lZi^ ^ 

Oil, ai^aN-i-2 

and ^^{{z}) = ({z})) otherwise. The q—difference operators = l,2,...,iV} are given by: 

( 2 . 11 ) = D*^[^\{zi,Z2,.:,ZN}-,{o^o,o:i,...,aN-ei,o:N-\-2}) 


aoq 


and 


( ^ (-l)l'^-|4>.({z})Af vr + 1 + 

*{z} = ' 0 *{^}{{zi,Z 2 ,.:,Zk}]{aQ,ai,...,ak+i,Zk+i}) for 


ctN-ei 


q^al 


k = l,2,...,iV- 1. 


Importantly, these operators essentially coincide with the ones given in [1108] . up to an overall factor and 
additionnal constant term chosen for further convenience. According to [11081 Proposition 4.5], the g—difference 
operators , fc = 1, 2,..., N are self-adjoint and mutually commuting: 

(2.12) [D*fj\D*®]=0 forall A, 1 € {1, 2,..., A}. 

Below, we will sometimes use an alternative formula for the g—difference operators introduced above as polynomials 
of {Ez^}. According to [11081 Proposition 4.2], by induction on N one obtains the following result: 


Proposition 2.1 (See [1108] . Proposition 4.2). Define Zg = ag and zat+i = aA-i- 2 - The N—variable q—difference 
operator can be written as: 


(2.13) 

where 

(2.14) 


j*(A) ^ QQg 

aN+i 






ag{q^ + 1 ) 


XgXN-\-l 


llfc=l Ok (z) 
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with the convention i/q = ua+i = 0. For j = 0,1, 


BY[z) = 1 + 


“i+i 


= (^1 - '^2+l^p+l.^ (^1 _ 


q^aj+iZj+iZj 


q'^a'j {q^ + 1)0;^ ’ ' V / \ 

V aj J \ ajZj J V aj J \ ajZj+iJ 

= Ij{Ij+i{B]’\z))), B-^’\z) = I,{By{z)), B^’-\z)=I,+^iB^’\z)) fork,I = 0,1, 

and 

b^Jz) = (1 - q^z^)il - q^zf), b]{z) = (1 - z^)(l - q^zf), bj^z) = WM). 


Remark 2.2. According to i2.11\) . • From Provosition \2. R it follows: 




^OtN+1 _ -I ^N+1 

^O^A ^+1 — 1 “I” 2 9 


- Y. c'.({z}). 

!/£{-!, 0 ,l}'^\{ 0 }'^ 


In view of the invariance of the normalized Gasper-Rahman polynomials (12.101) under the action of the involution 
f (see Theorem 12.21) . a ‘dual’ family of mutually commuting difference operators can be introduced. Let 

, g"’’,..., q--)[Etl,Etl,..., Etl] 

denote the associative algebra of difference operators with rational functions of g"^,..., as coefficients. 
Following [1108] . introduce the map b (which extends f) such that: 

u( \ u! \ CtoaN+lOiN+2q f ■ 1 o :\r I 1 u( ^ 

b(ao) = ao, b[aj) = - for j = 1, 2 ,..., + 1, b(aA+ 2 ) =-, 

OiN+2-j Oioq 

\j(^z ) = „2(wi+ra2 + ...+ww+i-,) 


aog 


-q 


h{E, 


= for j = l,2,...,N with = M 


The commutative subalgebra An of is generated by N algebraically independent difference operators {ID)|^^j|A: = 
1,2, ...,iV} which are defined as follows: 

Definition 2.5. Let k = 1,2,...,7V. The dual mutually commuting N-variable difference operators are defined by: 


(2.15) 


1 )W 


*(fe) 

1^1 


)■ 


By construction |II08| . the normalized multivariable polynomials of Gasper and Rahman (12.101) form a basis of 
the vector space Vx- Extending the well-known property of the Askey-Wilson polynomials, they also solve a family 
of bispectral problems associated with Az and An- 

Theorem 2.3 (See [1108] . Theorem 5.5). Let {a} € (C*)^"*"^ and k = 1,2, ...,A^. The normalized multivariable 
polynomial Q*^'^^({n}, {a;}, {a}) solves the following system of q—difference-difference bispectral problems: 

(2.16) ID*«QW(M,M,M) = (^g2^'= + ^g-2«'=)QW({n},{x},{a}), 

^ ^ V «o9 Ckfe+I J 

(2.17) I])[^^jQ(^)({n},{4,{a}) = (zw+i-fc + 2 )^+i_fc) (9^^^(W, {x}, {a}) 

where Rlk = ni + n 2 + ■ ■ ■ + Uk- 

For = 1, the above equations produce the well-known second-order g—difference equation (12.31) and three- 
term recurrence relations (12.41) of the Askey-Wilson polynomials, respectively. From the point of view of the 
representation theory, in the basis of Askey-Wilson polynomials, the operator defines a semi-infinite tridiagonal 
matrix. For N generic, in the basis of Gasper-Rahman polynomials, the corresponding matrix - which entries are 
determined by - enjoys a ‘block’ tridiagonal form, as we now show. 
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Lemma 2.1. The N—variable difference operator can be written as: 
(2.18) '''' " 


*S - E ( 


t^N -l—l'N •••L'2—L'a 1 —L' 2 ] JT'I^N T?l'N —l—l'N 77’L'2—L'S IT’l —1^2 
^0nin2---nN ^ni ^un-i 


I ^N—1 •••i^2 ^3 1 ^ 2 ] Tffl^N TT'l'NTpl'2—1^3 77’—1 —^^2 

-h ZCnin 2 ---nAr 

I -l—^N •••^2—^3 —^ 2 ] TTl'N Iff^N-l—^N J7'^2—^3 J7' — ^2\ 

-r ^CLnin 2 ---nM ^n 2 * ‘ ‘ ^un-i j ' 

Proof. Recall that the g—difference multivariable operator as defined by (|2.1ip can be alternatively written as 
f|2.13p . Now, observe 


b (0 = b(«. 


■eZ) = {EZEz\,m 


un riNj^iJ 

=Id 








_ J^l^l—1'2 T^i^2 — ^3 , , . fpl^N-l—l'N JT'l^N 


Decomposing according to the values vi = +1,0, —1, one obtains (12.181) . 


□ 


For practical applications, the explicit expressions of the coefficients entering in the decomposition ()2.18l) are 
needed. In the next subsection, they will be given explicitly for = 1 - in which case the coefficients entering in 
the three-term recurrence relation (12.41) are recovered -, as well as for N = 2. 

Remark 2.3. For N generic, the coefficients entering in i2.18\) are given by: 


(2.19) 

( 2 . 20 ) 
( 2 . 21 ) 

and 

[ 00 --- 00 ] 

^nin 2 -"nAr 


uWn L'AT-l—L-AT ■■■^2—^3 1 —*^ 2 ] _ 

^n\n2‘"nN — 


ai 


Jl^AT I^AT-l— ■■■1/2—1^3 —1 —^^ 2 ] 


^nin2 ■ ■ ■riN 


[l/N ■■■J/2—1^3 -^^ 2 ] 


^nin2-"nN 


‘^O'n+iO'N+2 

ai 

2aN+lOiN+2 

ai 

‘^o:n+io:n+2 


b {cA{z})U=+i ), 
b {CA{z})Z=-i) , 
b{CA{z})Z=o) 


ai 


‘2oiN+iaN+2 


+ 


o:n+io:n +2 


2 ai 


- E ( 


,[VN VN-I-VN ■■■V2 — V3 1-1/2] I [VN VN-l-fN ■■■V2—V3 -I-I/2] 
Cnin 2 ---RAc ' 


i-nin 2 -"niv 


i/e{-i,o,i}'v\{o}^ 


[l/jv 1/iV-l—1/iV ■■■1/2 —1/3 — 1 — 


^nin2---n],l 


1 - 1 / 2 ]^ 


Note that the expression of the last coefficient follows from Remark 12.21 

Now, according to Lemma [2d1 it is clear that doesn’t leave invariant the eigenspace of generated by 

{a:}, {a}). Precisely, the action of D|^j is as follows: 


Proposition 2.2. Let p € IN be fixed. Let Vp denote the eigenspace generated by the normalized mul¬ 
tivariable polynomials {Q*^'^^({'n}, {x}, {a})| ni + 112 + ... + = p}- On Vp, the N—variable difference operator 

acts as: 

(2.22) C Pp+i + Pp + Vp-i. 

Proof. For N generic, assume p = ni + 112 + ... + riAf is fixed. Observe: 

(2.23) b(X)Q^^H{n}, {x}, {a}) = + VN,n 2 + IZN-I - I^N, ...,nN + yi- 1^2, {a:}, {«}) G Pp+i/i- 

Using Lemma ETT] and the definition of Vp, the claim follows. □ 

Above arguments obviously extend to any operator with k = \,2,..., N —1. All these operators are mutually 
commuting semi-infinite block tridiagonal matrices. 
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2.3. New infinite dimensional modules for the q— Onsager algebra. We now endow the vector space 
C[xitX2, ■■■,xn] with a module structure of the g—Onsager algebra. A family of homomorphisms indexed by 
the integer /c = 1, 2,..., N can be exhibited as follows. 


2 


- 


(g 2 _ q-^f 


Proposition 2.3. The map defined by 
(2.24) Wq i-a X]\[-^.l—k, Wi 

is an homomorphism from Oq{sl 2 ) to 'Dz- 

Proof. First, consider the case k = N. Let denote the following linear transformation: 


(2.25) 

(2.26) 


W, 


W . 




w . 


C[a;i,X 2 , ...,XAr] I-A 

/ ^ 

/ ^ 


C[a:i,X2, ...,XAr] 
Xif, 

N, z 

"W P 


lTn)*W) ; 


To prove the claim, we begin with the first relation in ()2.1I1 . Let A* denote 


(2.27) 




(g2 - q-2)2 




Jg-M ' ^ u''o 

We now show A* = 0. According to the definition of the multivariable g—difference operator, recall that is 

a linear combination of with vi = 0, ±1. Then, observe: 

(2.28) {2x,f - (g2 + q-^)2x,Et"i2xi) + (Et,\2xi)f + (g^ - q-^f = 0, 

which implies A* = 0. 

We now turn to the second relation in (12.11) . Recall {Q(^)({n}, {x}, {a})}, DIat = 0,1,2,... form a basis of 
the vector space C[xi,X 2 , ■■■,xn]. With respect to this basis, according to Theorem 12.31 the operator is 

diagonalized by ({n}, {x}, {a}) with eigenvalues: 


(2.29) 


^*iN) ^ 1 / CtN+l 2mN Qog -2gtw^ 

2 V aog aN+i ) ' 


On the other hand, by construction, the action of on {Q^^^({n}, {x}, {a})} produces a linear combination 

of polynomials {x}, {a})} with OJIat = mi + ... + mN and DJIat = 0,1,2,.... Let denote the 

entries of the corresponding matrix. According to Proposition [521 it has vanishing entries for |9lAr — 97lAr| > 1 i.e. 
the operator acts as a block tridiagonal matrix in the basis {Q^^\{n}, {x}, {a})}- Let A denote the matrix 

representing 


(2.30) 


A = [Wr\ [Wr\ Wr^]g]g-i] + 


(g 2 - q-^Y 




with respect to the basis {Q^^\{n}, {x}, {a})}. We now show A = 0. To do this, we show that each entry of A 
is zero. According to the action of on the basis {Q^^\{n}, {x}, {a})}, the matrix A has {{n}, {m}) 

entry: 


T-C)’)"ww 


( 2 . 31 ) ((«;>';'>)» - (,»++ (ff'r+iP-fP) («; 

Recall = 0 for jOlAf — > 1- Observe that the first factor of the expression (12.311) is vanishing for 

|DlAf — 9^Ar| = 1 in view of the eigenvalues (12.291) . Obviously = 0 for |DlAf — WInI = 0. It follows 

A = 0. Both relations being satisfied, the claim is proven for k = N. For 1 < fc < A, using Theorem 12.31 similar 
arguments are applied from which the claim follows. □ 


Remark 2.4. The above arguments are identical to those in |Ter0321 p. 17-18] for N = 1. 
Let us now describe explicit examples of homomorphisms for N = 1,2. 
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Example 2.1 (See also [Ter032) . Section 5). For N = 1, on the vector space C[a;i] the elements Wq, \N\ of the 
q— Onsager algebra act, respectively, as = Xi and as the q—difference operator: 

(2.32) Wf) = ^ + 1 + , 

with 


$ 1 (^ 1 ) = 


(1 - - ^zi)(l - ^^ 1 ) 




Note that this operator coincides with the so-called Askey-Wilson second-order q— difference operator. 


Example 2.2. On the vector space C[a;i,a; 2 ] the elements Wq, W\ of the q—Onsager algebra act, respectively, as 
( 2 ) 

Wg ' = Xi and as the 2—variable q— difference operator: 

(2.33) = ^($ii(zi,Z2)A.iAz.+$ ll(^^^^2"')VzlV,, - $n(zi,Z2”')A.iV,,-$ii(zf\z2)V,,A,, 

+ ‘l’lo(2l! Z2)Az,, — $ 10(^1 22 )^zi + $01 (^1) Z2)Ai.Z2 ~ $0l(^lj ^2 ^)^22 + 1 + 

with 


0^3 \ 

q^al ) 


^ll{zi,Z2) 

^lo{zl,Z2) 

^m{zi,Z2) 


(1 - aizi)(l - f|zi)(l - f^2iZ2)(l - q^^,z^Z2){l - ^Z2){1 
" {l-zl){l-q^zl){l-zl){l-q^zl) 

(1 - aizi)(l - f|zi)(l - ^^Zi)(l - ^Zi) 
{l-zl){l-q^zl) 

(1 - a2Z2)(l - %Z2){1 - ^Z2){1 - ^Z2) 
(l-z|)(l-g2^|) 


:Z2) 


Note that using the definition (I2.11I1 . the explicit homomorphism (12.2411 for = 2 and fc = 1 is given by 
Wg^^ = X 2 and acts as (I2.32|) with as —>• Z 2 . 

For generic values of N, according to Lemma ETT] the spectral problem (12.171) gives a coupled system of three-term 
recurrence relations with respect to the integer (see Proposition l2.2l) . In the basis of normalized Gasper-Rahman 
multivariable polynomials (12.101) . in view of the homomorphism (j2.24l) for k = N, it implies that Wg can be written 
as a block tridiagonal matrix which entries are given by (I2.19L (I2.20p . (12.211) . Let us now describe some explicit 
examples for = 1, 2. 


Example 2.3. For N = 1, in the basis {Q^^\ni,xi,aQ,ai,a2,a3)\ni = 0,1,...} the operator Wq acts as a semi¬ 
infinite tridiagonal matrix denoted Wg^^. It coincides with the well-known operator that produces the three-term 
recurrence relations for the Askey-Wilson polynomials. Following the conventions above, on Q^^'^(ni,xi, {a}), 


(2.34) 

with 


W, 


( 1 ) 


acts as -f -f , 


,[o] 


^ni 


ai 


(1 - a2a3g'"0(l - 


2a2Q!3 


(l-^^4n,-2)(l-^^4n,) 


4 - 1 ] _ 0 ^ 20:3 




OL^a^ 


aL°l = 


2ai 

0 . 20.3 

2o\ 


(l_^^4n,-2)(l_^q4„,-4) 


0\ 


- foil] - c'-^l 
20203 ■ 
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Example 2.4. For N = 2, in the basis {Q^^\ni,n 2 ,xi,X 2 ,cio,ai,a 2 ,a 3 ,a 4 )\ni,n 2 = 0,1,...} the operator Wq 
acts as a semi-infinite block tridiagonal matrix denotedyV^\ On Q^^^{ni,n 2 ,xi,X 2 ,{a}), 


(2.35) 


with 


W, 


( 2 ) 


acts as 


5 ( 10 ] E + E 
_l_ E 

' ^n\n2 '^ni ' ^nin2 ^ri2 ' ^nin2 


'^nin 2 ■ 

^ ^711712 ^'^1^772 ^ ^ni7l2 ^71^ ^^2 ^ ^7li7l2 


-2 

712 


6 L'°i = 




fo[01] 

^711712 


.[- 10 ] _ 


.[ 1 - 2 ] _ 


Ql 


(1 - a3a49"^^)(l - - ^q4„i+2n2-2)(i _ ^^4„i+2n2)(l _ ^q2ni-2)(i _ ^^2ni) 


2a3a4 (1 _ ^54942-2)(1 _ ^g4942)(l - §94"i-2)(l _ ^g4ni) 

(1 - a3a4g''"=)(l - _ ^q2n2)(l _ ^^2„2+2 )(i _ 5-2ni)(l _ 


ai 


^Oq2-2n,) 


2a3a4 


(1 - ^g4m2-2)(l _ ^94042)(1 _ ^g2-4n,)(l _ ^g4-4n,) 


ai 


(1 - a3a49''^=)(l - 


_5[io] _;,[-i 2 ] 

^nin 2 ^711712 ’ 


2a3a4 (1 _ ^q4?n2-2)(i _ ^g4^n2) 

^304 (1 - - ^g4ni+2n2-2)(i _ ^^4m + 2„2 -4) _ ^2ni)(i _ gg 2 ni- 2 ) 

^ (1 - f|g401.-2)(l _ ^g4042-4)(l _ ^g4n.-2)(l _ ^g4n,-4) ’ 


0 : 30:4 


2q:i 


(1 _ ^54942-2)(1 _ ^g4>342-4)(l _ ^g2-4„i)(i _ ^g-4ni) 


n — ^3 „2<Xt2-2\/^ _ 03 „2942-2Wi _ 2942'\('1 _ £lfl2D42-2'| 

JO- 1 ] _ 03041^ 04^ ^ og*^ 


2q:i 


,[ 1 - 1 ] = 


,[- 11 ] = 


,[ 00 ] _ 


(l-^g4^n2-2)(l-^q49l2-4) 

“0 “0 

ai (1 - - a304g2"3)(i 


_ f.[- 10 ] _ Jl- 2 ] 

^711712 '•'711712 ’ 


2030:4 

0304 (l-g^-)(l-f9^--^)(l- 


(l-^g4ni-2)(l_^^4ni) 


_ ?,[ 10 ] _ ^[ 1 - 2 ] 

'^711712 '•'711712 ’ 


„ 2 ni —2 


Q;5Q;3a4 


)( 1 - 


„ 2 ni — 2 \ 


2 ai 

03 O 4 


(l_^^4n,-2)(l_^g4„,-4) 


_ r[-4°l - 51-42] 

'-'nin 2 '^711712 ’ 


ai 


_ „[ 1 - 1 ] _ „[-lll _ 5 [- 12 ] _ /,[ 10 ] _ ;,[ 01 ] _ [ 1 - 2 ] _ [- 10 ] _ [ 0 - 1 ] 

'^711772 '^7li7l2 '^7li772 '^77i772 '^7l\7l2 '^77,171-2 ‘^711712 '^71177,2 ' 


2 ai 20304 

More generally, realizations of the elements of the g—Onsager algebra in terms of g—difference operators are 
obtained in a straightforward manner using Proposition 12.31 Also, realizations in terms of difference operators are 
obtained from the following identification: on the polynomial j^;}^ |q;})^ 

(2.36) acts as ^®{n} ^^id acts as 

with p.l 8 [) . the expressions given in Remark 12.31 and p.29p . 

Note that for certain sets of parameters ag, oi,..., aAr 4 - 2 , the inhnite dimensional module of the g—Onsager 
algebra generated by the Gasper-Rahman polynomials may become reducible (see Section 3). Also, the conditions 
of irreducibility of this module are not studied in this paper. 

Combining previous results, we obtain one of the main results of this paper: 

Theorem 2.4. Let V denote the infinite dimensional vector space Vx- Let I/I/q, Wi he the standard generators of 
the q—Onsager algebra Ii2.1\) with q 1, where I/I/q, l/l/i act on V as i2.25\) . i2.26\) . Then, I/I/q, Wi induce on V a 
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module structure for the q— Onsager algebra with p = —{q^ — q ^)^/4. The Gasper-Raman multivariable polynomials 
{x}, {a})|{n} € (N)'^} form a basis of this module. 

In Section 3, we will consider subspaces of parameters ...,aM +2 for which the module generated by the 

Gasper-Rahman polynomials becomes reducible. In this case, the semi-infinite dimensional matrices associated 
with truncate to finite dimensional ones. 


2.4. Raising and lowering operators. In the context of orthogonal polynomials with one or several variables, 
the determination of the raising and lowering operators is an important issue. For the multivariable Macdonald 
polynomials associated with a root system of type A„_i, see for instance [KN96] . For the Askey-Wilson polynomials, 
examples of raising and lowering operators have been constructed either starting from the g—difference-difference 
bispectral equations or using the one-variable version of the double affine Hecke algebra |Sa07] . In this subsection, 
by analogy with [Sa07] . we extend the construction to the multivariable Gasper-Rahman polynomials (I2.10|) in 
a straighforward manner. Let p € IN he fixed. Let Vp denote the subspace of C[xi,X 2 ,..., xat] generated by the 
normalized multivariable polynomials {^}i {o^})| ni -I- n 2 -I- ... -I- njv = p}, as in Proposition 12.21 Below, 

±(fe) 

we introduce N distinct pairs of difference operators which act as: 

(2.37) C for all fc = 1, 2,..., N. 

Note that raising and lowering operators for the Askey-Wilson polynomials are recovered by setting A^ = fc = 1 in 
the formulae below. 


Proposition 2.4. LetWQfWi be defined as Define for a given set of integers {n} = {ni,n 2 , ...ynk): 


(2.38) = 




_ jn+ei] 

n k-i (n*{i) _ o*{i) \ fa*h) _ n*{h \ I \ ( 

/=! ^{n+ei})\^{n} ^{n-ei})/ \r{n±efc} ^{n^ek}J {n±ek} ^{n} ) 


q*(0 n*{l) 


*{k) 


g*(fc) 




n*{k)\ 


w, 


(fc) 


For any N > k and any {n'} such that n' = Ui for I < i < k, one has 


(2.39) 

(2.40) 








:LO:.:„o-iiF;-iQW(K},{x},{a}). 


Proof. Let us first prove Equation (12.391) . In the definition (12.381) of define F as the result of the action of 

the term on {x}, {a}), which is determined by ()2.36p with (j2.18() for N ^ k. By construction, 

through the action of (resp. on F, all the terms for which Uk = 1 (resp. Vk = —1) 

vanish. If Vk = 0, a factor of the form 0^^^^ — (resp. 9*^^^ — appears in front of the corresponding 

terms. Therefore, the I = 1 term in the numerator of ()2.38l) projects out all the terms in F for which i//; = ±1 and 
leaves the terms for which = 0 unchanged up to a factor 

in front of each of the remaining terms. Then, looking at the action of the operators with Z = 2 in the numerator 
of (j2.38|) . which have the same structure as those for Z = 1, on the remaining Vk = 0 terms of F, we similarly find 
that it projects out the terms for which Vk-i — Vk = Vk-i = ±1 and leaves the terms for which Vk-i = 0 unchanged 
up to a factor. An immediate induction on 1 < Z < Zc — 1 then shows similarly that only the terms for which t'i = 0 
for 2 < i < k are not projected out, and those are left unchanged up to a factor. Therefore there now only remains 
three terms from F: and acting on {x}, {a}). Finally, the second part of 

(12.381) acts on the remaining expression as follows: (wj^^ — 0*^l^y) / projects out the term in 

whereas (w}^^ — projects out the term in , and again both leave the term 
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in Ern. unchanged. Thus, dividing the final expression by 

fc-i 

1=1 

in the end, only one term has survived from T. It is given in the r.h.s of (12.391) . The proof of the other equation 
(12.4011 is similar. □ 

Let us mention that in the context of quantum integrable systems, the solution of the ‘inverse problem’ - i.e. the 
explicit relation between local operators and the elements of the non-Abelian algebra that ensures the integrability 
of the model - plays a central role in the analysis of correlation functions, for instance. The explicit construction of 
the raising and lowering operators (12.381) in terms of the fundamental elements of the g—Onsager algebra provides 
such a solution. It should find applications in the analysis of models generated from the g—Onsager algebra, for 
instance the open XXZ spin chain. 

3. BiSPECTRALITY AND THE RELATION WITH TRIDIAGONAL PAIRS 

As mentionned in the introduction, for the g—Askey scheme of orthogonal polynomials in one variable the 
bispectral property has a natural interpretation within the framework of Leonard’s theorem |BT84) and Leonard 
pairs |Ter87l[Ter031) . Namely, for a pair of elements A, A* that act on an irreducible finite dimensional module, are 
diagonalizable, and satisfy the so-called Askey-Wilson relations, it is known that there exists two basis with respect 
to which the two matrices representing A, A* are diagonal (resp. tridiagonal) and tridiagonal (resp. diagonal). In 
such examples, the overlap coefficients between the two basis are given by the Askey-Wilson polynomials evaluated 
on a discrete support. First examples appeared, for instance, in |Z91j . 

Tridiagonal pairs are generalizations of Leonard pairs [Ter04) , in which case one assumes that degeneracies may 
occur in the spectrum of the two matrices under consideration. Provided the elements A, A* satisfy a pair of 
relations called the tridiagonal relations (see (j3.3l) . (13.4|) below), it implies the following: in the basis in which the 
first matrix is diagonal with degeneracies, the other matrix takes a block tridiagonal structure. Furthermore, there 
exists another basis with respect to which the first matrix transforms into a block tridiagonal matrix, whereas 
the second one transforms into a diagonal matrix with degerenacies. In the literature, it is now well-understood 
that the representation theory of tridiagonal algebras - in particular the g— Onsager algebra jni) - is intimately 
connected with the theory of tridiagonal pairs [Tern32) . as we will recall below. For instance, finite dimensional 
irreducible modules of the g— Onsager algebra for which A and A* are diagonalizable have been classified in |IT09] 
for g not a root of unity. As an application, tridiagonal pairs of g— Racah type (see ‘case F below) over C have 
been classified for g not a root of unity. For an algebraically closed field and no restrictions on g, note that a 
classification of tridiagonal pairs is given in |INT10) (see also |IT081) 1. However, to our knowledge, the connection 
with the theory of orthogonal polynomials has remained, in general, an open problem. The purpose of this Section 
is to show that the bispectral property of the Gasper-Rahman polynomials finds a natural interpretation in the 
theory of tridiagonal pairs. From that point of view, given a TD pair, the overlap coefficients between the two basis 
are given by the normalized Gasper-Rahman polynomials (12.101) evaluated on a discrete supporiF^. In this context, 
a correspondence between the finite dimensional tensor product evaluation modules of the g— Onsager algebra 
considered in |Bas06l IIT09| and the factorized structure of Gasper-Rahman polynomials (j2.5|) arises naturally. 

3.1. Tridiagonal pairs and tridiagonal algebras. In the present subsection, we recall the results taken from 
[Ter93[ llTTer991 [Ter032[ INTOQj . Let K be an arbitrary field unless otherwise noted. Let U denote a vector space 
over K with finite positive dimension. For a linear transformation A : U ^ U and a subspace W C U, we call W 
an eigenspace of A whenever W ^ 0 and there exists d £ IK such that W = {u & U \ Au = 6 u]. 0 is the eigenvalue 
of A associated with W. A is diagonalizable whenever U is spanned by the eigenspaces of A. 

Definition 3.1 (See [ITTer99] . Definition 2.1). Let U denote a vector space over IK with finite positive dimension. 
By a tridiagonal pair (or TD pair) on U we mean an ordered pair of linear transformations A : U ^ U and 
A* : U ^ U that satisfy the following four conditions. 

(i) Each of A, A* is diagonalizable on U. 


^^For applications to quantum integrable systems generated by the g—Onsager algebra, scalar products of the Hamiltonian’s 
eigenstates will be expressed in terms of the overlap coefficients. 
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(ii) There exists an ordering {Up}p^Q of the eigenspaces of A such that 

(3.1) A*Up C Up-i + Up + Up+i 0<p<d, 

where C/_i = 0 and Ud+i = 0 . 

(iii) There exists an ordering {U*}p^Q of the eigenspaces of A* such that 

(3.2) Au; c u;_^ + u; + o<p<s, 

where = 0 and = 0 . 

(iv) There does not exist a subspace W of U such that AW C W, A*W C W, W 0, W U. 

We say the pair A, A* is over K. We call U the underlying vector space. 

Note 1 . According to a common notational convention, for a linear transformation A the conjugate-transpose of 
A is denoted A*. We are not using this convention here. In a TD pair A, A*, the linear transformations A and A* 
are arbitrary subject to (i)-(iv) above. 

Let A, A* denote a TD pair on U, as in Definition 13.II By [ITTer99l Lemma 4.5] the integers d and 5 from (ii), 
(iii) are equal and called the diameter of the pair. An ordering of the eigenspaces of A (resp. A*) is said to be 
standard whenever it satisfies (13.111 (resp. (13.21) 1. Let {Up}p^Q (resp. {D*}p^q) denote a standard ordering of the 
eigenspaces of A (resp. A*). For 0 < p < d, let Op (resp. 9*) denote the eigenvalue of A (resp. A*) associated with 
Up (resp. [/*). According to |Ter032[ Theorem 4.4], for K = C the eigenvalues Op and 9*, p = 0,1,..., d, can take 
three different forms; 

• Case 1: the g—Racah type (with b, c, b*,c* 0) 

Qp = a\- bq^P-‘^ + cq’^-^P, 
e*p = a* + b*q^P-‘^ + c*q‘^-^P-, 

• Case 11: the Racah type (with c, c* ^ 0) 

Op = abpcp{p — l)/2, 

0 *p = a*+ b*p + c*p{p - l)/ 2 ; 

• Case 111: the Bannai-lto case 

0p = a + b{-l)P + cp{-l)P, 

0* = a*+b*i-l)P + c*pi-l)P. 

Here a, a* ,b, b*, c, c* are scalars and q ^ 0, cf^ f^l, q'^ ^ —1. In this case, we say that A, A* is a tridiagonal pair 
of g—Racah (case I) or Racah (case 11) type |NT09l Theorem 5.3]. Note that the cases listed above are the ‘most 
general’ cases in which the main free parameters are nonzero. Other examples for which some of the parameters 
are set to zero can be found in [Ter033) . 

The theory of tridiagonal pairs is closely related with the representation theory of tridiagonal algebras. Tridi¬ 
agonal algebras are defined as follows: 

Definition 3.2 fSee [Ter032| . Definition 3.91. Let (3,^,^*, p, p* denote scalars in M.. WedefineT = , p, p*) 

as the associative f^-algebra with unit generated by two elements A, A* subject to the relations 

(3.3) [a,A^A* - ^AA*A-^A*A^--/{AA*-^A*A)-pA*] = 0, 

(3.4) [a*,A*^A-^A*AA*+ AA*^--f*{A*A-^AA*)-p*A] = 0. 

We refer to T as the tridiagonal algebra (or TD algebra) over K with parameters /3, 7 , 7 *, p, p*. A and A* are called 
the standard generators ofT. 

In the literature, note that the relations (I3.3L(I3.4I) are called the tridiagonal (TD) relations. In case I, if 
(3 = q^ q~‘^, 7 = 7 * = 0, p = p* ^ 0, the parameter sequence is said to be ‘reduced’. In this case, the TD relations 
coincide with the defining relations of the g—Onsager algebra m with, for instance, Wq —>■ A, Wi —> A*. 


Remark 3.1. For ,0 = 2, 7 = 7 * = 0, p = p* = 16 the TD relations coincide with the Dolan-Grady relations 
[DoG82| . For (3 = q^ q~‘^, 7 = 7 * = 0, p = p* = 0, the TD relations coincide with the q—Serre relations of 
Uqish). 
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It is known that every TD pair conies from an irreducible finite dimensional module of a tridiagonal algebra. 
The following theorems establish an explicit connection between the representation theory of tridiagonal algebras 
(finite dimensional irreducible modules for q not a root of unity) and the theory of tridiagonal pairs, as defined in 
Definition 13.11 


Theorem 3.1 (See [Ter032) . Theorem 3.7). Let A, A* denote a TD pair over K. Then there exists a sequence 
of scalars p* taken from IK such that A, A* satisfy the tridiagonal relations VS.Al . i3.4\ )- The parameter 

sequence / 3 , 7 , 7 *,p, p* is uniquely determined by the pair if the diameter is at least 3. 


Example 3.1 (See [Ter032) . Lemma 4.5). Let A, A* be a TD pair of q—Racah type (case I). Suppose a = a* = 0 
and be = b*c*. Then, A, A* satisfy the defining relations of the q—Onsager algebra i2.1\} with Wq —^ A, \N\ —>■ A*, 
p = —bc{q^ — q~^Y. 

Theorem 3.2 (See [Ter032] . Theorem 3.10). Let (3,'y,j*,p,p* denote scalars in K, and assume q is not a root of 
unity, where q^ + q~‘^ = j3. Let T denote the TD algebra over K with parameters fd, 7 , 7 *, p, p* and standard gener¬ 
ators A, A*. Let U denote an irreducible finite dimensional T-module and assume each of A, A* is diagonalizable 
on U. Then A, A* act on Lf as a TD pair. 


According to some observations in |Ter032[ p. 11], let A, A* denote a TD pair on U with parameter sequence 
/3,7,7*,P,P*. Let / denote the identity and assume (3^2. Then, the TD pair 


(3.5) A = p^ ' (i + 7(/3-2)-i/) , A*=p^^(^p* + 


T 


*2 


2-/1 


(i* + 7*(/3-2)-i/) 


has a reduced parameter sequence: the TD pair (13.51) satisfy the defining relations of the q—Onsager algebra (the 
so-called g—Dolan-Grady relations) ( 12 . 11 ) with Wq —>■ A, Wi —>■ A*. It is easy to show that such transformation is 
invertible. As a consequence, provided /3 ^ 2 , finite dimensional irreducible modules of the tridiagonal algebra with 
7 ,p (resp. 7 ,p) both not 0 and of the g—Onsager algebra are isomorphic. For this reason, in the next subsection 
we simply focus on the construction of finite dimensional modules. The special case q = 1 for which /3 = 2 will be 
treated separately in Section 5. 


3.2. Finite dimensional modules. In view of the homomorphism (12.241) ioi k = N together with (I2.36|) . a com¬ 
parison between the structure of the spectrum of a tridiagonal pair of q—Racah type for q ^ 1 (called ‘Case F above) 
and the bispectral problem (I2.16L (12.171) strongly suggests that finite dimensional modules of the q—Onsager alge¬ 
bra m can be derived by studying certain quotients of the infinite dimensional modules based on the normalized 
Gasper-Rahman multivariable polynomials (12.101) . considered in the previous Section. 

Considering bispectral relations of the form (12.3L (12.41) . it is well-known that finite dimensional modules of 
dimension 2j -\- 1 based on the Askey-Wilson polynomials (12.21) can be easily constructed provided the parameters 
a,b,c,d are chosen such that the coefficienTi 62 ^ = 0, and that 2 : is restricted to a discrete support. In this case, 
the q—difference equation (12.31) becomes a three-term recurrence relation in addition to (EH), and both recurrence 
relations terminate. Then, the corresponding bispectral problem is interpreted as the spectral problem associated 
with a Leonard pair. Such a technique can easily be extended in order to construct finite dimensional modules 
based on the normalized Gasper-Rahman polynomials (I2.10L in view of the factorized structure (|2.5I) . To this end, 
restrictions on the parameters ai ,..., aAr+i and a discrete support of the variables zi, 22 , zn are considered. 

Let Wo,Wi be the standard generators of the q—Onsager algebra m with q ^ 1. Recall that Wo,Wi as 
defined in Proposition [231 induce a module structure on C[a;i, 0 : 2 ,..., xat]. To construct finite dimensional modules 
of the q—Onsager algebra, we proceed as follows. Let hi, fi 2 ,..., un be positive integers. Assume the A^-variables 
are now restricted to the following values: 

(3.6) _l_ ^2 -I - + hk for fc = 1, 2,..., iV. 

Oik 

Then, for the following discussion we introduce the notation: 

(3.7) aW({n}, {n}, {a}) = QW({n}, W, 2 a' 


12 


'By construction, note that the coefficient cq = 0. 
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According to the results of the previous Section, on the restricted support generated by (EH), one induces 
operators having the same action on {h}, {a}) as xi and the same action on ({n}, {x}, {a}) as 

For simplicity, we will not change their names and identify them again with Wq, Wi (see Proposition 13.II 
below). They are still diagonalizable. The (degenerate) eigenvalues of Wq (resp.Wi) are, respectively, given by: 


(3.8) 




( “1 „-2TtAr 1 

“Ar+iaA+2 „2^nA,^ 

a*iN) 1 1 

^ “09 - 2 ' 71 n I 

“JV+I 2<nN\ 

2 ^ 

1 Q 1 

\0'N-\-l0'N-\-2 

y 1 5 

Oil / 

‘'{"I 2 ' 

VQiiV+i 

^ 1 
o^oq ) 


and the corresponding eigenfunctions are £j(^(*({n}, {h}, {a}). 

Then, the vector space is spanned by the eigenspaces of either Wq or Wi. Let p G N be fixed. With respect to 
Wo,Wi, two different families of ordered eigenspaces can be constructed, respectively, from 

= Span{n(^(’({n}, {h}, {a})|ni + ... + uat = p}, {n} fixed, 

or = Span{0(^(’({n}, {h}, {a})|ni + ... + riAT = p}, {n} fixed. 

Up to now, note that the set of integer variables {n}, {h} has not been restricted. Now, provided the parameters 
Qfi,..., ciAr+i satisfy certain conditions, finite dimensional modules can be constructed as follows. 


Proposition 3.1. Let ji,j 2 , be positive half-integers or integers and define dN = 2(ji+j2 + ---+jAf)- Assume 

there exists parameters oi, 02 ,..., oat+i such that: 

(3.9) = k = l,2,...,N. 

Otk 


Then, the maps defined by 

(3.10) Wo^xi, 


IVl 


*{N) 


(g2-q-2)2 


yield a family of homomorphisms from Oq{sl 2 ) to any of the vector spaces ^ fixed, and 

also to any Vp^'^ with n fixed. These are finite dimensional modules of the q— Onsager algebra 112.1\) . 


1 ^ f ./\^ ^ 

Note that, as stated at the top of this page, although Xi and are actually distinct operators for each of 

the target vector spaces, we leave the same notation for all of them since they are induced from the same operator. 


Proof. First, we have to show that is a finite dimensional module of the g—Onsager algebra (12.11) . According 

to the results of the previous Section, on the space generated by (13.71) the action of the standard generators Wq, Wi 
is characterized by (12.361) . Observe that the element Wi acts on ^pZo'^p^^^ a diagonal matrix with entries 
such that 2flN = P- Given the ordering of the eigenspaces {Vp^^^}p=o consider the action of 

the element Wq, using (12.361) with (12.181) . Observe that Wq acts as (13.2F To prove that the module is finite 

dimensional, it is sufficient to show that the coefficients (I2.19L(I2.20|) . (I2.21I) are such that: 

0 if vn = 1 and ni = 2ji, 

VN-k — VN+i-k = ^ and nk+i=‘^jk+i^ 

IZ 2 = 0 and nN = 2jN, V 2 =-1 and tiat = 2jAr, 2jAr - 1, 
VN-k — VN+i-k = 2 and Uk+i = 2jfe+i, 2jfe+i — 1, 

0 if i^Af = 1 and ni = 2ji, 

VN-k — VN+i-k = 1 and Uk+i = 2jfe+i, 

VN-k — VN+i-k = 2 and Uk+i = ‘^jk+i^‘^jk +1 — 1, 

0 if i^Af = 1 and ni = 2ji, 

VN-k — VN+i-k = ^ and nk+i=‘^jk+i^ 

V 2 = —1 and uai = “^jN, 

VN-k — VN+i-k = 2 and Ufc+i = 2jk+i^‘^jk+i — 1, 


/o 1 1'1 uWn VM-l-VM - l-V 2 \ 

Onin2---n;v 


[VN VN-l-VN ■■■V2-V3 -1-1^2] 
Gnin2 • • -niv 


[vN VN-l-VN ■■■V2-V3 -V 2 \ _ 


^n\n2---nN 
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for all A: = 1,2,iV — 2. To show that, one uses the formulae given in the Appendix: given the conditions on the 
ratios of the parameters (j3.9|) . one finds: 

b(S^’°(z)) = 0 and b{Bl~^{z)) = 0 for UN+i-k =‘ 2 jN+i-k, k = l,...,N-l, 

= 0 for m = 2 ji, 

b{Bl~^{z)) = 0 for riAT+i-fc = 2jAr+i_fc, 2jAr+i_fc - 1, k = l,...,N-l. 

As a consequence, using (j2.14p the coefficients ()2.19p . (j2.20p . (|2.21ll automatically satisfy (13.lip . 

Secondly, we have to show that is also a finite dimensional module of the g—Onsager algebra Recall 

(13.61) . On one hand, according to (12.241) for k = N, the element Wq acts on as a diagonal matrix with 

entries On the other hand, introduce the ‘dual’ difference operator: 


(3.12) 


_ I 

l{ao->-ao; «>=-> ° "^^1' 


Then, with respect to the indices {h} it follows that Wi act on according to (13.ip . Now, using the explicit 
expressions of the coefficients (12.141) . one finds that the coefficients of ‘ 6 —type’ characterizing (13.121) are vanishing 
for fik = jN+i-k with fc = 1,2, As a consequence, Wq, Wi act on as a diagonal matrix and an block 

tridiagonal matrix, respectively. The claim follows. □ 


Remark 3 . 2 . Note that 


(3.13) 


dim 



dim 


/ dM N 

0v;w 


Kp—O 


N 


l[i2jk + l). 

k=l 


Example 3 . 2 . For N = 1 , let j be a positive half-integer or integer. Define 

= Span{Q^^^(ni, hi, ao,«!, 0:2, 0:3)}, ni fixed, 
= Span{Q^^^(ni,hi,ao,Q:i,Q:2,a3)}, ni fixed, 

with 


(3.14) 


ao = 


•-B-l 


ai = q a2 = q 0:3 = 5 


where /3,/3* are arbitrary scalars. The finite dimensional module (resp. of the q— Onsager algebra has 

dimension 2j + 1. The eigenvalues of Wq, l/l/i, are distinct and given by: 

gi}) = i ^^/3+2i-2«i ^-p-2j+2hiS^ ^ ^*(1) ^ 1 ^^/3*+2j-2ni ^-/3*-2i+2ni ^ ^ 

Note that using (I3.14p in (12.341) . one can check that 6211 = 0 and cq = 0. 


Example 3 . 3 . For N = 2 , let ji,j2 be positive half-integers or integers. Define 

= Span{ 0 (^^(ni,n2,hi,h2,00,01,02,03,04)1711 +712 =p}, ni,n2 fixed, 

= Span{£l^^^(77i, 712 , hi, h2, oo, oi, 02, 03, 04)|hi + h2 = p}, ni,n2 fixed. 


with 


(3.15) 


oo = q^ 


01 = 5 


02 = 5 ^ 


03 = 5 t 


04 = 5 


where 13 , 13 * are arbitrary scalars. The finite dimensional module (resp. of the q—Onsager algebra 

has dimension (2ji + l)(2j2 + 1). The 2(ji + J 2 ) + 1 degenerate eigenvalues of Wq, W\, are given by 

= - (nP+'^(n+32)--^3l2 , -p-2(ji+j2)+2<h2\ n*i2) ^ 1 / fl*+2(ji+j2)-29r2 , -N-2Ui+32)+2<n2\ 

^hih2 2 ^ ^ / ^ 1^2 2 ^ ^ / ’ 
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Note that using (I3.15|) in (12.3511 . for any ni,n 2 one can check that: 


„[-10] 

^0 tI2 

= 0, 

JO-1] 

^niO 

= 0, 

Jl-2] . 

^niO ■ 

= 0, 


c 

^[10] 

"2jin2 

= 0, 

^[Ol] 

Oni2j2 

= 0, 

51 - 12 ] 

0™i2j2 

= 0, 

‘LJt. = 

0 , 

“niO 

= 0, 

[-11] 

« 0 n 2 

= 0, 

“2iin2 

= 0, 

AS=o- 



[1-2] ^ Q 

2jin2 

%n2 ~ 


Example 3.4. For N generic, let jk = 1/2 for all k G {1,2,..., JV}. Define 

y*{N) ^ Span{0(^)({n},{n},{a})|91jv =ri, fit fixed, 

V^^'> = Span{5^^^({n},{h},{a})|0lAr =p}, Uk fixed, k = 1,2,...,N, 

with 

(3.16) = ai=g“^, ak+i=q~^~^ for k=l,...,N, aN+ 2 =q~^, 

where /?,/?* are arbitrary scalars. The finite dimensional module (resp. ) of the q—Onsager algebra has 

dimension 2^. The N + 1 degenerate eigenvalues of Wq, l/l/i, are given by 

nW _ 1 f B+N-2<hN , „-/3-IV+2aijv^ n*{N) _ 1 / /3*+Ar-2miv , „-/3*-IV+2miv\ 

^{h} -2V ® J’ {"} “ 2 ® / ■ 

Based on the previous results, explicit examples of tridiagonal pairs can be constructed as follows. Let U 
denote a vector space over C with finite dimension. Let : U ^ U and : U ^ U denote two linear 

transformations depending on parameters ao,ai, ...,aN +2 satisfying Proposition 13.II Assume U has a basis 

(3.17) {V'(ni,n 2 , ...,n7v)|0 < Uk < jk,k = 1,2, 
with respect to which, on i/(ni, n 2 ,..., u-at), 

(3.18) acts as ^®{n} acts as 

Assume U has a ‘dual’ basis 

(3.19) {ip{ni,n2, ...,hAr)|0 < fiN+i-k < jk: k = l,2, ...,N} 
with respect to which, on (p{fii,fi 2 , ...,n 7 v), 

(3.20) acts as and acts as 

The conditions (/), {ii) and {in) of Definition 13.11 being fulfilled, assume U is irreducible. Then, is a 

tridiagonal paii[{j of q—Racah type, as defined in Definition 13.11 By analogy with the Askey-Wilson polynomials 
that are interpreted as the overlap coefficients between dual basis for Leonard pairs |Z91[ITer^7llTerO^ . it follows: 


Theorem 3.3. Let U be the finite dimensional module of the q— Onsager algebra a which admits the two ‘dual’ 
basis ^3.17^ , i3.19\) . Assume U is irreducible. Then, Wq, l/l/i aet on U as a TD pair of q—Raeah type defined as 
The multivariable polynomials defined on a discrete support {lj*^^)({n}, {h}, {a})} are the overlap 
coefficients between the two basis. They solve the difference-difference bispectral problem: 

Wr0W({n},{h},{o}) = 0WdW({n},{h},{a}), 

wW5W({n},{n},{a}) = 0*^S^dn}, {h}, {a}), 

where are given by \3.8\) . 


the condition (zn) is relaxed, then VVq^\ is called a mocfc tridiagonal pair. For details, see | IT082| . 
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In the case = 1 of one-variable polynomials, one recovers the usual duality for Askey-Wilson polynomials 
found in Leonard’s theorem [BT84) . 

As mentionned at the beginning of this Section, finite dimensional irreducible modules of the g—Onsager algebra 
(EH) were classified in |IT09[ Theorem 1.15], using an embedding of the g—Onsager algebra into the Uq{£{sl 2 ))-^oop 
algebr [1T091 Proposition 1.13]. In the analysis, tensor product evaluation modules of the g—Onsager algebra 
[1T091 Section 1.4] (see also [BasOBl Lemma 2.5]) 


(3.21) 


V{2jN,VN) 0 • • • 0 V{2ji,vi) 


were considered, where Vj € C are usually called the evaluation parameters. From a general point of view, for g 
not a root of unity a TD pair of g—Racah type is afforded via the embedding of the g—Onsager algebra into the 
Uq{£{sl 2 ))-loop algebra by (13.211) . It has diameter d = 2ji -|- ... -|- 2jN and the dimension of the vector space on 
which Wo,Wi act follows from [IT091 Proposition 1.24] with A —>■ l,£i —?► 2A,n —?► N, which result coincide^ 
with (13.131) . Note the factorized structure of the Gasper-Rahman polynomials (12.51) by analogy with the factorized 
structure of finite dimensional modules of the g—Onsager algebra considered in |IT09[ lBas06| . 

To conclude, let us mention that finite dimensional modules of the form (13.211) have found applications in 
the solution of the open XXZ spin chain with generic integrable boundary conditions |BK07| . For the family of 
tridiagonal pairs discussed in [BasOB] . up to a normalization, the eigenvectors in each basis have been determined 
explicitly. In particular, they possess a factorized structure (see Proposition 3.3 and 3.7 of [BasOB] L parallel to the 
factorized structure (|2.5I) . 


4. The g—D olan-Grady hierarchy and spectral problem revisited 

As mentioned in the Introduction, there is a rather large class of quantum integrable models on the continuum 
or lattice whose local integrals of motion can be written in terms of the elements of the so-called g—Dolan-Grady 
hierarchy HH), an Abelian subalgebra of the g—Onsager algebra (12.11) . For the integrable models that fall in this 
class, finding the spectrum and eigenstates of the Hamiltonian dm relies on studying the spectral problem of 
the elements {l 2 fc+i|fc = 1,2,...} using the representation theory of the g—Onsager algebra. For instance, in the 
literature, the two-dimensional Ising and superintegrable Potts models have been studied in details for g = 1 using 
the explicit relation between the Onsager algebra and a fixed-point subalgebra of SI 2 (under the action of a certain 
automorphism of s /2 |Dav 90j). For g yf 1, the explicit relation between the g—Onsager algebra and a certain coideal 
subalgebra of Uq{sl 2 ) has been used to analyze the finite open XXZ spin chain |BK07| or its thermodynamic limit 
analog for various types of boundary conditions |BB12l IBKol4j . According to the size of the system - finite or 
semi-infinite -, finite or infinite dimensional representations (g—vertex operators) of the g—Onsager algebra have 
been considered. 

In the previous Sections, modules of the g—Onsager algebra - either infinite or finite dimensional - have been 
constructed in terms of Gasper-Rahman multivariable polynomials. The purpose of this Section is to reformulate the 
spectral problem for the mutually commuting quantities {bfe-i-i} as a system of g—difference or difference equations 
(the case g = 1 will be considered separately in Section 5). It implies that eigenfunctions of Hamiltonians of the 
form EH) can be expressed as combinations of Gasper-Rahman multivariable polynomials. As an application, a 
g—hypergeometric formulation of the eigenfunctions of the open XXZ chain with generic boundary conditions is 
given for generic parameters, see ((131) . Also, a new derivation of Nepomechie’s relations is proposed and generalized, 
see (14.61) . 

4.1. Infinite dimensional modules. If we assume that the spectrum of li is non-degeneratcF^. common eigen¬ 
states of bfc-i-i are uniquely determined from solving the spectral problem of li. For generic parameters uio,u!i,g±, 
we assume such hypothesis holds true. Then, the action of this operator with respect to the basis of Gasper-Rahman 


^"^Note that a different embedding can be considered, see | Baa06l Proposition 2.2]. 

^^We especially thank Paul Terwilliger for pointing out Proposition 1.24 and other results of [ITQ9| to us. 
^^It is the case for the open XXZ spin chain, see |BK07| . 
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polynomials is derived as follows. Denote: 


<j)0 


E 

E 


CA{z})U=±iEZ---eZ, 


Ty /r IM i ^OiN+1 

E v\\z})\i,-^^—qE ■ ■ ■ E ^ + xoXjv_i_i. 


Using the results of previous Sections and the definition of Gi,Gi in terms of Wo,Wi given in (11.31) . on the 
multivariable polynomial basis, the elements Wq, Wi of the g—Onsager algebra (EU and their hrst ‘descendants’ 
Gi,Gi act, respectively, as: 


(4.1) 



= 



aoq 

, 

2aA+i 


aoq 

((g3 - g-i)(zi$+:B,, + + (g - q-^){z, + , 

AaN+i 

or 

aoq 

((g - g-3)(^-ici>+A,, + zi<l>-u;/) + (g - q-^)iz, + . 

4aAr+i 


As a consequence, on the multivariable polynomial basis the first element li = ujq^Nq + wiWi +g+Gi +g-Gi of 
the g—Dolan-Grady hierarchy (EH) acts as the multivariable g—difference operator: 

= 2^ " 2 " ^ (g+g^i + g-g-^u 

^("^1 + i(g+ + g-)(g - g-')(zi + ‘J>°- 


Alternatively, in the multivariable polynomial basis, the above operators EE) can be written as semi-infinite 
block tridiagonal matrices, whose entries are explicitly expressed in terms of (I2.19L (12.201) and (12.211) . By straight¬ 
forward calculations combining all expressions, on the multivariable polynomial {{n}, {x}, {a}). 


r(^) 


acts as 


{ 1 ^ 2 ,!^3,. 




VM-l—VN 
i^nin2 • • -n-iv 


•1 — 1^2] rpl'N JT'l^N-l — l'N . . . Tp^ — ^2 


1 — - 1 — 1 ^ 2 ] J^i^N fP^N-l — ^N J^l^2—l^3 fp — ^ — ^2 

"T ^mn2---nN ^n2 ‘ ‘ ' ^un-i ^riN 

_|_ A\y^ ^N-l—l^N - ^“ 2 ] -pl'N — fpl^2—i-'3Zp — l^2\ 

“T ^nin2---nN ^n2 ' ’ ’-^riN-i -^riN 


where 


(A 0\ ^N-i—^N •••1 — 1 ^ 2 ] 

y^.Z) *^nin2---nN 


^Wn ^N-1—^N -1 —^^ 2 ] 

^rL\n2‘‘‘ni^ 


A[i^N ^N-l—l^N -^^ 2 ] 

>^nin2---niv 


I (g^ ~ g ( Qfpg 2inw-i I „ ofAf+i 2inw+A\.k;v 
I CJq “t \d— Q ' 9-\- 0. I J^nin2-"nAr 

V 2 \ qat+i aoq // 

Lo + +5+^g-^"'"+^))cLrn;”4 

+ (<^0 + (g++ g-)(g — g ^ 


For generic parameters g±, the above coefficients are non-vanishing so that there is no multivariable 

polynomial subspace left invariant under the action of the semi-infinite block tridiagonal matrix representing li. 
Corresponding eigenfunctions of the spectral problem: 

jWvl/W({x}, {a}) = aWv|/W({x}, {a}) 


(4.3) 























BISPECTRAL <j-HYPERGEOMETRIC STRUCTURE AND q - OA’S INTEGRABLE MODELS 


23 


are, in general, not polynomials. However, given the multivariable polynomial infinite dimensional basis of the 
vector space, they can be written as: 

(4.4) vI/W({4,{a})= ^ 

{n}^]N^ 

where the coefficients /|„} satisfy a coupled system of three-term difference equations, with respect to the integer 
For generic parameters and q not a root of unity, the above eigenfunctions diagonalize the full g—Dolan-Grady 
hierarchy. As a consequence of (11.11) . the Hamiltonian is diagonalized by (14.41) too. 

Example 4.1. For N = 1, recall {a} = {oo, <ai, 0 : 2 , 0 : 3 }. The eigenfunctions of the Hamiltonian H.l|) are given 
by {a}) = {o:}) where the coefficients fm satisfy the three-term recurrence relation 

(f-i=0,fo = l): 

/ni-iSi\Li + fn,+iClf[+i + fm (A™ - Ai) = 0 for ni = 0,1,... 

Example 4.2. For N = 2, recall {a} = {ao,ai,a 2 ,cx 3 ,a 4 }. The eigenfunctions of the Hamiltonian are 

given by ^^=0 ^ 2 ,a^i, 0 : 2 , {a}) where the coefficients fm,n 2 satisfy the coupled 

recurrence relations (f-i,n 2 = /m.-i = /ni+ 1,-2 = 0/' 

f f fplOi] I f 

Jni-l,n2‘-’ni-l,n2 J "i ."2, 712 -1 .''il + l."2-2‘'’ni + l,712-2 

I f ^[—10] . n ^[0—1] . r ^[1 — 2] 

“•“7^1 + 1,n2^ni + l,712 + 712 + 1 ^“l)^2+2'-^7ij^ —l,n2+2 

+ /ni-l,n 2 + l'4lLj_p„2 + l + /ni + 1,n2-l'4lLi + l,712-1 F /’ll ,"2 ('4lni°,n2 ~ ® /'^^ 111,712 = 0,1,... 

Note that according to Theorem 12.11 provided the parameters {a} satisfy the conditions (12.81) . scalar product 
between the eigenstates (14.41) can be computed in a straightforward manner and expressed in terms of g—shifted 
factorials (see (E3). 


4.2. Finite dimensional modules. Let V be defined as in Proposition 13.11 For generic parameters LUo,uJi,g±, 
using the bispectral property, an eigenstate of li (and more generally of H given by ()l.ll) l can be written in terms 
of Gasper-Rahman polynomials defined on a discrete support, either 

dN dN 

(4.5) T(^)({fi},{a}) = {+,{+) or 4'^ ({n}, {a}) = ^ /{ii|n(^)({n}, {h}, {a}), 

= l T!n = 1 


where the coefficients /{„} (resp. f^ay) satisfy a coupled system of three-term difference equations, with respect 
to the integer 91iv (resp. Ofjv). The explicit form of these equations follows from the action of li on the bispectral 
polynomial basis: the system of recurrence relations determining fy^} is associated with the coefficients ()4.2I) . 
whereas the ‘dual’ system of recurrence relations determining fyjiy is associated with: 


*[!7jv VN-l—VN "T —!72] 


®7li7l2---7liV 
r,*[uN I^N-l — VN 

■'h-ih2---nw 


-1-772] 


= Wl 


= Wl 


(g^ -q ^) /„ aAf+iaAf-l-2 „29fw-l-l 


[9+ 


(g2 - 9-2) 


ai 

ai 


+ 9 -- 


ai 


o:n-\-io:n-\-2 


-g 




7T,l7l2---72iV 


-*^ 2 ] 

'^hih2---nN 


(^ 

ujo9{a}S^,o + (wi -b {g+ + g-){q - g“^)fi'{ft})a~+ 


ci;Ar-i-iaA-i-2 


^-2 T!n-I-1 ^ Q^A'-|-lQ:A-|- 2 ^ 2 ^w-l^^^*[’ 7 iV un-i-un 

~ ai 


-'rii7l2"-7lAr 


UN UN-1—I'M - ^ 2 ] 

h2---nN 


Here, the coefficients ^ follow from the transformation (13.121) in (12.181) . 

Remark 4.1. Let V be the finite dimensional module of the q—Onsager algebra, of dimension 2^, defined as in 
example \3.4\ -Any of the states |.).5[ ) provide an eigenstate of the Hamiltonian of the open XXZ chain with generic 
boundary conditions. 


1 ^ 2 ] 

7 

— U 2 
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4 . 3 . Invariant subspaces and Nepomechie’s relations. For special relations among the parameters wq, wi, g±, 
multivariable polynomial subspaces that are left invariant under the action of the element li can be constructed. 
Let V be defined as in Proposition 13 . II In particular, max(iXlAr) = max(9lAr) = d^. Let P,P* be integers. Define 


W+ = Span{0^'^^({n}, {h}, {a})|0 < tllAr < {«} fixed, 

W- = Span{0^^^({n}, {h}, {a})|0 < OTat < P*}, {n} fixed, 

W+ = Span{Q(^^({n}, {h}, {a})|P + 1 < fTlAT < max(fltv)}, {h} fixed, 
W_ = Spcin0^\{n},{h},{a})\P*+ l<mN < max(iiljv)}, {n} fixed. 


According to the structure of the coefficients of the block tridiagonal matrices associated with b, it follows: 


c w+ 

if 



c W- 

if 



i^^w+ c W+ 

if 

.. , (9^-9-^) 
W 0 + 2 

{g- 

c W- 

if 

.. , (9^-9-^) 
UJ1+ 2 

( 5 + 


OtN+l 
aN+iaN+2 2P*+i 


9+ - 9 

aoq 


■)=»■ 


ai 


+ 9-- 


ai 


o:n+io:n+2 


= 0 , 


aoq 


+ 9+ - 9~ 

0 ! N +1 


= 0 , 


Q:i _2P*_i Q!Ar+lQ!Ar+2 2P* + 1 

- 9 +9- - 9 

Q!Ar+iQ;Ar+2 cti 


= 0 . 


On each of the invariant subspaces, eigenstates of the Hamiltonian (ED can be systematically constructed as 
truncations of any of the two possible combinations (14.5|) . For explicit examples, see [BK071 Section 3.2]. 

In the literature on quantum integrable systems with boundaries, special cases of the above relations among 
the parameters in (14.61) first appeared in the algebraic Bethe ansatz analysis of the open XXZ spin chain with 
non-diagonal boundary conditions. Sometimes refered as ‘Nepomechie’s relations’ |N03[ eq. (1.4)], they usually 
characterize the family of constraints between the left and right integrable boundary conditions of the system for 
which the standard Bethe ansatz approach can be applied. Within the q—Onsager approach of the open XXZ 
spin chain, they were independantly identified from a representation theory perspective in [BK07] . Importantly, 
it provided a rigorous proof of completeness of the spectrum of the Hamiltonian. Note that the above relations 
also occur in other approaches, see for instance [FaKN13] . as well as they characterize a class of solutions to the 
boundary g—Knizhnik-Zamolodchikov equations |SV13j . Here, we would like to stress that the relations (14.6p 
generalize the results of [BK07] : they characterize the existence of a class of invariant subspaces of the Abelian 
subalgebra of the g—Onsager algebra. As a consequence, they should also arise for integrable models associated 
with higher spin representations, for instance. 


5. The special case q = I 

In this Section, new infinite and finite dimensional modules of the Onsager algebra [044] with defining relations 
(12.11) for q = 1 are constructed in terms of the multivariable Krawtchouk polynomials introduced by Tratnik |T89] 
(see also [GerI07| L The analysis is extended in a straightforward manner to more general tridiagonal algebras with 
defining relations (13.3|) . (13.41) for /3 = 2. Restricting the parameters allows to derive explicit examples of tridiagonal 
pairs of Racah types (Case H), as briefly discussed. 

First, we recall the duality and bispectrality of the multivariable Krawtchouk polynomials [T891 [GerI07] . Al¬ 
though the notations may slightly differ, the basic material is taken from |GerI07] to which the reader is refered 
for more details. In terms of the Gaussian hypergeometric function, the one-variable Krawtchouk polynomial is 
defined as: 


kn{x,a,b) = (- 6 )„ 2 Fi 


—n, —X 1 
-b ’a 


n G N and {x, a, b) G 


with {y)n = T{y + n)/r{y) a function which can be well defined for all y real and n integer. Following |T891IGerI07j . 
A"—variable generalizations of the Krawtchouk polynomial can be introduced. 
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Definition 5.1 (See |T891IGerI07| '). Let oi ,aN and M he nonzero real parameters. The normalized N—variable 
generalization of the Krawtchouk polynomial is defined by: 

1 Af 

(5.1) k^^'>{{n},{x},{a};M) =- - Y[knfixj;- 

where the notation 


M — ‘TIn + — Xj-i), 


JCj — Xi “t" X2 


a: j j — 0 ^ 


1-21,-1 

21, = 01 + Oi2 + • • • + Q!,, 2to = 0 


is used. 


By analogy with the g—deformed case, an inner product on the space of polynomials V can be introduced. For 
f^g^V, it takes the form: 

(5-2) (/,£/)= f{x)g{x)p{x). 

With respect to this inner product, for a positive integer oq and An < 1 the multivariable Krawtchouk polynomials 
(EH) are orthogonal on the domain JF = {{a:} £ N-^IXat < M} for the weight p{x) given by: 

= {n-xnv. n ^ ( 1 ^^) ■ 

Remark 5 . 1 . The normalized multivariable Meixner polynomials are obtained from fOD through the substitution: 

ak = - — - and M = — s . 

Cl + • • • + cn — 1 


Following [GerlOTj . multivariable Krawtchouk polynomials are known to diagonalize two different commutative 
algebras of difference operators in the variables {x} and {n}, respectively. Let 

denote the associative algebra of difference operators with rational functions of Xi,X2, ...,xn as coefficients. A 
mutually commuting family of difference operators which are diagonalized by the multivariable Krawtchouk poly¬ 
nomials (EU has been constructed in [GerlOTj . Here, they are denoted = 1,..., N}. 


Definition 5.2 (See [GerlOTj . Section 5.4). The N-variable Krawtchouk difference operators = 1,2, ...,N} 

are defined as: 


(5.3) 

where 


dl 




{x} 






, XN, 


ClN+l-l aN-f-2-l 


OiN 


1 — 2lAr_i 1 — ^N-l ' 


1 — 2lAr_z 


; M — Xn-1 


N N 

CN{xi,X2,...,XN,ai,a2,....,aN;M) = {Xn - M)Y. ajExj + (SIat - 1) x^Ef^ 

i=l i=l 


N 

- Y, ^i^j^XjEf} - Y. + ^n{M - Xn). 

l<i^j<N i=l 


The algebra generated by these difference operators can be understood as a limiting case of the commutative 
algebra Az of the g—difference operators considered in Section 2. According to [GerlOTj . the difference operators 
^*{x} are mutually commuting: 

(5.4) [dl*{xldl*{x}]=0 for all Z,jG{l,2,...,iV}. 

As explained in [GerlOTj . the normalized multivariable Krawtchouk polynomials (15.ip are invariant under the 
following change of variables 

“Ar+ 1 -, (1 - An) ,, ,, 

Xj riAr+i-,-, rij XN+I-,-, a, -!> (i _ ■)(l-AAr ■)’ M ^ M. 























26 


PASCAL BASEILHAC AND XAVIER MARTIN 


In particular, it implies the duality relation [GerlOTl Theorem 5.6]: 

(5.5) W, {«}; M) = {i}, {a}; M) 

By analogy with the case g 1, an isomorphism bi between the associative algebras of difference operators in 
the variable {a:} and {n} that extends the above map can be introduced: 


bi(M) = M, bi(a,) = 


OiN+l-j{^ — -4 at ) 


for j = 1,2,..., Af, 


bi(a;j) = UN+i-j, bi{Ea:^) = for j = 1,2, 

Then, a ‘dual’ commutative algebra generated by a family of difference operators {d[^^^y\l = 1, ...,iV}, follows from 
(15.31) through the action of bi. According to the definitions of [GerlOTj . one has: 

Definition 5.3 (See [GerI07| l. Let Z = 1, 2,..., N. The dual mutually commuting N-variable Krawtchouk difference 
operators are defined by: 


(5.6) 


<} = &!(<})• 


Example 5.1. For k = N, the dual N-variable Krawtchouk difference operator takes the form: 

N N 

(5.7) dl[^l = {UlN-M)Y,HaN+i-i)En,+bi{AN-l)Y.n,E-^ 

N 


- ^ bi(aAf+i-,)?tj-A„,A„/ + giAf + bi(.4jv)(M - 9djv) - bi(aAf+i-0?t,. 


i=l 


As a consequence of the duality relation, it follows |GerI07] : 

Theorem 5.1 (See [GerI07| . Theorem 5.7). Let ({a}, M) € and I = 1,2, ..., N. The normalized multivariable 

polynomial K^^\{n},{x},{a};M) solves the following system of difference-difference bispectral problems: 

(5.8) d7*®i?W({n},{4,{a};M) = (fyiN - mN-i)K<^^\{n},{x},{a}-, M), 

(5.9) dlfl^k(^H{n},{x},{ay,M) = XN+i-ik^^H{n}, {x}, {a}; M). 

For A^ = 1, the equations (EH), dSH) produce the well-known difference equation and three-term recur¬ 
rence relation of the one-variable Krawtchouk polynomials [KS96] . respectively. In particular, in the basis gen¬ 
erated by the Krawtchouk polynomials, the operator defines a semi-infinite tridiagonal matrix. For N 

generic, in the basis of multivariable Krawtchouk polynomials (EID the structure of the matrix associated with 
dl^^y is determined as follows. Let p G N be fixed. Consider the subspace Vp generated by the polynomials 
{K^^'>{{n},{x},{a}-,M)\ OIat = p}, and act on it with (15.7F The action of the operators A„., E~^ and ErnE~^ 
with i j in (15.71) on Vp sends it to Ip+i, Vp-i and Vp, respectively. As a consequence, 

(5.10) dl[^lVp C Kp+i + Kp + Kp_i. 

For N generic, we conclude that the matrix associated with is semi-infinite block tridiagonal. For A^ = 1,2, 
note that the explicit expressions of for A^ = 1 (/ = 1) and A^ = 2 (/ = 1,2) can be obtained from [GerI071 
Appendix A]. Above arguments obviously extend to any operator with I < N. 


5.1. Infinite and finite dimensional modules of the Onsager algebra. The defining relations of the Onsager 
algebra are given by (ED with (7 = 1. Let us endow the vector space K[a;i,X 2 , ...,xn] with a module structure of 
the Onsager algebra. By analogy with the analysis in Section 2, a family of homomorphisms indexed by the integer 
Z = 1, 2,..., N can be exhibited as follows. 
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Proposition 5.1. The map defined by 

(5.11) I/I/q i-A XAr_|_i_;, l/l/i i-A p I—>■ 1, 

is an homomorphism from Oi{sl 2 ) to T)'^. 

Proof. To show the claim for I = N, one essentially follows the same steps as in Proposition 12.31 The proof uses 
the fact that 

(5.12) xj - 2xiEt,\xi) + -1 = 0, 

(DIa, ± if - 2 {^n ± 1)01n + 01^ - 1 = 0, 

together with the observation that the matrix representating Wg is semi-infinite block tridiagonal (it acts as (Ih.lOp l 
in the basis which diagonalizes Wi. Similar arguments are used to show the claim for 1 < I < N. □ 


It follows: 

Theorem 5.2. Let V denote the infinite dimensional vector space R[a;i, 12 ,..., xjv]. Let H/g, IVi be the stan¬ 
dard generators of the Onsager algebra 12.11) with q = I, where Wq, Wi act on V as S5.111 . Then, 1/1/g, Wi 
induce on V a module structure for the Onsager algebra with p = 1. The multivariable Krawtchouk polynomials 
K^^\{n},{x},{a};M)\{n} G (N)^;M G M} form a basis of this module. 

Let us finally construct finite dimensional modules of the Onsager algebra in terms of multivariable Krawtchouk 
polynomials. On one hand, according to the expression (I5.3|) . if Xj = 0 observe that the contributions from the 

downward shift operator Ef^ disappears. On the other hand, let J G N. Suppose Xi = J for a given i. Then, it is 

possible to cancel all contributions associated with the upward shift operators in the first and third terms of (j5.3p 
as follows: we require both that Xn < J and M = J. Indeed, with these two conditions, note in particular that 
when Xi = J, Xn = J and all the other Xj with j i are necessarily 0 which implies that the first and third term 
of (lOl) disappear. Because of the simple duality relation between Xi, E^, and Ui, En ,, these conditions clearly also 
restrict the Ui and in the same way. It follows: 

Proposition 5.2. Let J be a positive integer, and assume M = J. Define 

Vj^'^ = SpanlKT^-^^dn}, {a;}, {a}; J)|XAr < J} for IHn ^ J fixed, 

(5.13) = Span{^W({^}, < J} for Xn < J fixed. 

Then are finite dimensional modules of the Onsager algebra. 

Note that the corresponding tridiagonal pair has a reduced parameter sequence: the matrices associated with 
Wg, Wi satisfy the defining relations of the Onsager algebra (i.e. (EH) with q = 1, p = 1). The TD pair is of Racah 
type (Case II), and the spectra have a simple structure with a = a* = c = c* = 0 and b = b* = 1. 


5.2. Generalization to tridiagonal algebras with fi = 2. By analogy with the analysis of the previous Sections, 
it is straightforward to show that the multivariable Racah polynomials of Tratnik [T89| generate infinite dimensional 
modules of the tridiagonal algebras with defining relations (13.311 . (13.41) for the special case (3 = 2. Recall that the 
one-variable Racah polynomials are defined aQ 


r„(a;; a, b, c, d) = {a 1)„(6 -I- d -I- l)„(c -I- 1)„ 4 F 3 


—n, n-I-a-|- 6 -I-1, —a;, a;-I-c-I-d-|-1 
a-|-l, 6 -|-d-|-l,c-|-l ’ 


Let Coj Cij •■•I Cn-i -1 be real parameters. Let N,M be positive integers. The normalized multivariable Racah poly¬ 
nomials are given by |T89] (see also [GerI07) h 

RN{{n},{x},{Q,M) 


(5.14) 

where 


R^({n},{a:},{C},M) = 


- Co)<yiN nf-i(Cfc+l - (k) 


N 


RN{{n},{x},{C},M) 


^rik { 1 1 “1“ Cfc Co 1 5 Cfc+1 Cfc 1 1 

1 1 Cfc ^fc+l) 




the literature, the notation rn(A(ai); a, b, c, d) with A(ai) = x{x + c + d + 1) is sometimes used |KS96| . 
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with VIn < M. With respect to the inner product (15.2|) . the polynomials i?Ar({n}, {x}, {(}, M) are orthogonal on 

(5.15) F = {{x} G N^|0 < xi < X 2 < • • • < xjY < M}. 

The explicit expression of the weight p(x) follows from [GerI071 eq. (3.8)]. In the parametrization above, it reads 
(xo = 0,XAr+i = M): 


1 '^', tt r(Cfc+i - Cfc + 2;fe+i - Xfe)r(Cfc+i + Xfe+i + Xfe)-p-|- , r,.„ ^ 

= n —(x,„-xo'r(a+i+x.«+xo— n(c> + 


fc =0 


fc=l 


In |GerI07) . using the duality property of the multivariable polynomials (15.141) . two commutative families of 
difference operators have been introduced. They are diagonalized by Tratnik’s multivariable Racah polynomials. 
Using the notations in [GerI07] . let us denote: 

(5.16) = and + Cn+i-i) for l = l,2,...,N, 

where 2f and £” are defined, respectively, in |GerI071 eqs. (3.24),(2.6) or (3.4)] and |GerI07l eq. (4.11)]. By 
analogy with Theorem 12.31 it was shown: 

Theorem 5.3 (See [GerI07| . Theorem 4.6). Let ({C}, M) £ and 1 = 1,2,..., N. The normalized multivariable 

polynomial i?Ar({n}, {x}, {C}, M) solve the following system of difference-difference bispectral problems: 

(5.17) d7;i'Ji7^(M,{x},{C},M) = A*Wi7^(M,{x},{C},M), 

(5.18) cff|'i|i7^(M,{x},{C},M) = Ag}i?iv(W,{4,{C},M). 
where 

(5.19) A^i'j = -(0+i-Co)^i-^i(9ti-l), 

(5.20) A^2} = —(1 + CAI+l-O^^Af+l-/ ~ ~ !)• 


Infinite dimensional modules of tridiagonal algebras can be constructed as follows. We endow the vector space 
K[xi, X 2 , ..., xjv] with a module structure of a tridiagonal algebra T(2, 7 , 7 *, p, p*), as in Definition 13.21 A family of 
homomorphisms indexed by the integer 1 = 1,2,..., N can be exhibited as follows. 

Proposition 5.3. The map defined by 

(5.21) A a|2p 7 “2, 7 * i-A — 2, 

P Ca+i-; ~ P* ^ (Co — 0+1 + 1)^ ~ 1. 

is an homomorphism from T(2, 7 , 7 *, p, p*) to T)'^. 

Proof. To show the claim for 1 = N, one essentially follows the same steps as in Proposition 12.31 The main 
difference here is the structure of the two-variable polynomials associated with the l.h.s of (13.3|) . (13.41) . respectively. 
Similarly to (I2.28|) . (12.311) . these polynomials are vanishing due to the structure of the spectra (|5.19l) . (15.201) . 
provided the parameters 7 , 7 *,p,p* are fixed to the appropriate values. Similar arguments are used to show the 
claim for 1 < 1 < A^. □ 

It follows: 

Theorem 5.4. Let V denote the infinite dimensional vector space M[xi, X 2 , ..., xa]- Let A, A* be the standard 
generators of the tridiagonal algebra T(2, —2, —2, p, p*), where A, A* act on V as \5.21\) . Then, A, A* induce on V 
a module structure for T{2, — 2 , — 2 , C^_|_i_; — 1, (Co ~ Ci+i + 1)^ ~ !)• The multivariable Racah polynomials 
i?(-^)({n}, {x}, {C}, M)|{n} £ (N)^;M £ R} form a basis of this module. 

To conclude this Section, let us observe that finite dimensional modules of the tridiagonal algebra T{2, 7 , 7 *, p, p*) 
based on the multivariable Racah polynomials (15.141) can be easily derived by considering certain quotients of the 
space of parameters {(},M. To this end, it is sufficient to identify the set of constraints on the parameters - by 
analogy with Proposition 13.11 - using the explicit expressions for and £” taken from |GerI07[ eq. (3.24)] and 































BISPECTRAL <j-HYPERGEOMETRIC STRUCTURE AND q - OA’S INTEGRABLE MODELS 


29 


[GerI071 eq. (4.11)], respectively. The analysis being straighforward, we skip the details. Note that in the basis 
generated by the corresponding multivariable Racah polynomials (I5.14p defined on a discrete support, the elements 
A, A* act as a tridiagonal pair of Racah type with spectra of the form (15.191) . (15.201) . This is in perfect agreement 
with the expected structure of the spectra that follows from the theory of tridiagonal pairs of Racah type (Case 
II). 


6. Perspectives 

Besides its elegance, for the class of quantum integrable models generated from the g—Onsager algebra, to have 
a g—hypergeometric formulation should be of practical importance. It should give a fresh look on existing solutions 
following from different frameworks such as the algebraic Bethe ansatz or separation of variables approaches but 
not only. We think it could pave the way to circumvent or reconsider some of the difficulties inherent to the 
structure of already existing analytical solutions. Below, we list a conjecture and some open problems. 

Irreducibility criteria: We conjecture that the modules built on Gasper-Rahman polynomials are irreducible for 
generic parameters ao,ai,aN +2 and q, meaning that they do not take specific values such as those of Proposition 

[Q 

Open XXZ spin chain: For g ^ 1, one of the main motivation for constructing a g—hypergeometric basis 
associated with the g—Onsager algebra originates in the problem of solving quantum integrable systems with generic 
integrable boundary conditions, the XXZ open spin chain to cite the simplest example. On one hand, recall that 
for certain regimes of boundary parameters or generic parameters, through the Bethe ansatz |N031 ICLSW021IYZ071 
IFiKlOj , its modified versions |Bel41 IBeP141 IGrl41 IZLGYSW14) , functional alternatives [Ga071 IGal41 IGYSWl^ 
or the separation of variables approach [FaKNlSl IKMN14) . eigenstates and eigenvalues are expressed in terms 
of (Bethe) roots of highly transcendental Bethe equations. It appears that an explicit characterization of the 
Hamiltonian’s eigenfunctions as polynomials offers the possibility of studying the correspondence between Bethe 
roots and solutions of algebraic eouationJ^. On the other hand, in the thermodynamic limit, either from the 
algebraic Bethe ansatz setting |KKMNST07| or the g—vertex operator approach |JKKKMW94l[BKoI4| . correlation 
functions and form factors are always given in terms of multiple integral representations. Although formally 
satisfying, these are difficult objects to handle for practical applications. From the point of view developed here, 
the g—hypergeometric formulation (14.51) of the eigenfunctions opens the way to reconsidering scalar products 
of eigenstates in terms of Gasper-Rahman multivariable polynomials |GR041| . Whether corresponding integral 
representations could be related with the expressions derived within the g—vertex operator approach or algebraic 
Bethe ansatz is an interesting question. 

Higher rank q—Onsager algebras and orthogonal polynomials. Higher rank generalizations of the g—Onsager 
algebra (EH have been introduced in the literature [BB091 Definition 2.1]. They can be realized in terms of certain 
coideal subalgebras of Uq{'g) for any Kac-Moody algebra g |BB091 Proposition 2.1] (see also 1X121 1. For the simplest 
example g = 1 and g = such a generalization is called the Onsager algebra, first introduced in |UgIv95| . 

In view of the results presented here, an interesting problem would be to identify the family of multivariable 
orthogonal polynomials that would be higher rank generalizations of those of Tratnik (g = I) or Gasper-Rahman 
(g 7 ^ 1). Indeed, corresponding infinite and finite dimensional modules of the generalized g—Onsager algebras will 
provide an explicit g—hypergeometric basis for a class of quantum integrable models, among which one finds open 
spin chains with higher rank symmetries. 

Non-polynomial solutions and q— Virasoro: As an extension of Bochner’s original theorem, the bispectral problem 
for tridiagonal doubly infinite matrices has been considered in the literature. In this case, solutions are intimately 
related to the Virasoro algebra but are no longer necessarily polynomials [GH961] . In the context of g—special 
functions and the g—Virasoro algebra, a similar extension has been considered too [GH9621IGH97] . In view of 
these, it is clear that an extension of the bispectral problem ()2.I6|) . (|2.17|) for V > 1 can be considered similarly, 
see for instance |GerIII| . These works may establish a relation between the g—Onsager algebra and an algebraic 
structure extending the g—Virasoro algebra. 

Some of these problems will be considered elsewhere. 


^®Recall that zeroes of Askey-Wilson polynomials satisfy Bethe equations associated with a ‘I’-site open XXZ chain, as pointed out 
in IWZ95I . 
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APPENDIX: Useful formulae for the coefficients 


• For j = — 1: 



(1 “^+2-fcg4^n„. 

“o 

_fc+2nAr+i_fc —2^ 


/ “Af+2-fc 

O' 

_fc+2n^^4-^ -fc —2 W-| 

a2 9 

Oq 

a; 

KBl'-\z)) = 

2 2 
/-I ^Y+2 —fc ^2n^/4.i_fc—2^/i ^0 

V-^ 2 ^ 

^N+l-k 

/ V-^ 2 y 

^N+l-k 

II 
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(1 - q-‘ZnN+l-U^^l^l 
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2 y 
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KBl'-\z)) = 

2 2 
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V" q2 
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J\- 2 
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II 
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For j = 0,N: 
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